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Abstract 

We classify weak-scale extensions of the Standard Model which automatically 
preserve its accidental and approximate symmetry structure at the renormalizable 
level and which are hence invisible to low-energy indirect probes. By requiring the 
consistency of the effective held theory up to scales of Aeg ~ 10^® GeV and after 
applying cosmological constraints, we arrive at a hnite set of possibilities that we 
analyze in detail. One of the most striking signatures of this framework is the 
presence of new charged and/or colored states which can be efficiently produced in 
high-energy particle colliders and which are stable on the scale of detectors. 
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1 Introduction 


The electroweak (EW) hierarchy puzzle suggests that new physics (NP) degrees of freedom 
should appear around or not much above the EW scale. Hence the search for NP is a clear 
target being vigorously pursued by the LHC experiments. On the other hand, indirect searches 
for NP using flavor and CP violating observables have already probed NP scales up to 10® GeV 
(c.f. [1-3] for recent reviews). Thus, the overall excellent agreement with the CKM paradigm 
predictions suggests a large mass gap above the EW scale. Perhaps even more strikingly, 
searches for baryon (B) and lepton (L) number violating processes at low energies suggest that 
these accidental quantum numbers of the standard model (SM) are good symmetries of nature 
up to scales of the order of 10^® GeV (c.f. [2]). 

Explicit models of TeV scale NP need to resolve the apparent conflict between these two 
sets of expectations by postulating exact or approximate symmetries which in term forbid or 
sufficiently suppress the most dangerous contributions to flavor changing neutral currents, CP 
violation, as well as B and L changing processes. These include explicit B, L or their anomaly 
free combinations, discrete space-time symmetries including C, P, CP but also new internal 
symmetries like R-parity (or R-symmetry) in supersymmetric (SUSY) extensions of the SM, 
KK-parity in extra-dimensional setups, as well as abelian or non-abelian horizontal (flavor) 
symmetries. 

At the heart of these problems is the fact that generically, extending the SM particle content 
will either (i) break some of the SM accidental symmetries, and/or (ii) introduce new sources of 
breaking of the approximate SM symmetries, which in general will not be aligned with existing 
SM symmetry breaking directions. Examples of the first kind include B and L. Flavor, CP and 
custodial symmetry of the Higgs potential fall into the second category. 

Consider the SM as the renormalizable part of an effective field theory (EFT) 

£ = 411*’ + E . (1) 

d>4 

where only SM fields appear as dynamical degrees of freedom in £, and d denotes the canonical 
operator dimension. Assuming 0{1) coefficients in the EFT operator expansion, currently all 
experimental evidence in particle physics can be accommodated by such a generic theory with 
a very large cut-off Agg ~ 10^^ GeV.^ This particular scale is intriguing since it can account 
for both the observed neutrino masses suggesting the presence of L violating as well as 
null results of all flavor, GP and B violation probes constraining One may thus ask 

the following well defined question. Which extensions of the SM particle content with masses 
close to the EW scale (i) form consistent EFTs with a cut-off scale as high as 10^® GeV, (ii) 
automatically preserve the accidental and approximate symmetry structure of the SM and thus 
do not require the introduction of additional protective mechanisms in order to remain viable 
LHC targets in light of negative search results of the numerous indirect probes, and (in) are 
cosmologically viable? 

In the present paper we explore such possibilities by adding to the SM EW (and possibly 
color) multiplets, requiring that their SM gauge quantum numbers alone forbid all renormal¬ 
izable interactions which would break any of the SM approximate or accidental symmetries. 

^Cosmological observations suggest that most of the mass in the observable Universe cannot be accounted for 
by known forms of matter. The possibility of particle dark matter within our setup is briefly discussed below. 
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In particular, in the SM the global Qp = U(3)^ flavor symmetry of quarks and leptons is only 
broken by their respective Yukawas (and the gauging of hypercharge). In order not to introduce 
new sources of breaking, one should only consider Qp singlet operator extensions of the SM. 
In Sect. 2 we list all d < 3 operators involving quark and lepton fields that transform nontriv- 
ially under Qp and demand that the new degrees of freedom do not couple to any of these at 
the renormalizable level. Since both B and L are subgroups of the above prescription also 
automatically preserves these accidental symmetries of the SM. Furthermore, in most cases this 
also ensures the absence of new sources of breaking for both custodial and CP symmetries at 
the renormalizable level. The exceptions where new breaking can arise in the scalar potential 
are discussed in Sect. 2.2.1. 

The relevant Lagrangian is restricted only by imposing the SM gauge and Lorentz invariance, 
not by new symmetries. Finally, such a theory is assumed to represent a consistent description 
of nature up to the cut-off scale Aeff. In particular, we require that all of the marginal couplings 
(in particular the SM gauge couplings) remain perturbative up to Aeg. The physical idea 
behind this requirement is that a Landau pole might be associated with the emergence of some 
new, generic dynamics that will break the accidental symmetries of the SM at scales above the 
Landau pole. As we show in Sect. 2.4, this condition (together with the cosmological constraints 
on stable charged particles) limits the size of the new representations and leads to a finite list 
of possible SM extensions. 


Spin 

A 

Qlp 

C^decay 

dim(C>decay) 

K~ZZ[Gey] 

0 

(1,1,0) 

0 

XHW 

3 

> mpi {gi) 

0 

(1,3,0)^ 

0,1 

XHW 

3 

> mpi (^fi) 

0 

(l,4,l/2)t 

-1,0,1,2 

xHWW 

4 

> mpi (gi) 

0 

(l,4,3/2)t 

0,1,2,3 

XH^H^H^ 

4 

> mpi (gi) 

0 

(1,2,3/2) 

1,2 

xHHe, yti/tg^eA Df^x^ea^e’^ 

5 

> mpi (gi) 

0 

(1,2,5/2) 

2,3 


5 

> mpi (gi) 

0 

(1,5,0) 

0,1,2 

XHHH^H\ ylE^^lE^^, x^H^H 

5 

> mpi {gi) 

0 

(1,5,1) 

-1,0,1,2,3 

X^HHHW, yyyti/ti/t 

5 

> mpi {gi) 

0 

(1,5,2) 

0,1,2,3,4 

X^HHHH 

5 

3.5 X 10^® (c/i) 

0 

(1,7,0)* 

0,1,2,3 

X^H^H 

5 

1.4 X 10^® (c/ 2 ) 

1/2 

(1.4,1/2) 

-1 


5 

8.1 X 10^® (c/ 2 ) 

1/2 

(1,4,3/2) 

0 

xiH^H^ 

5 

2.7 X IQis (c/i) 

1/2 

(1,5,0) 

0 

xiHHH\ 

6 

8.3 X 10^^ (^ 2 ) 


Table 1: List of new weak-scale uncolored states y which can couple to SM fields at the 
renormalizable level without breaking Qp, and which are compatible with cosmology and an 
EFT cut-off scale of Aeff ~ 10^^ GeV. The possible electromagnetic charges of the LP in the 
multiplet are denoted by Qpp, while (Pdecay denotes the lowest dimensional operators responsible 
for the decay of y. States with Y = 0 are understood to be real. In the last column, the Landau 
pole has been estimated at two loops by integrating in the new multiplet at the scale of the Z 
boson mass m^, while the symbol in the bracket stands for the gauge coupling, 5 'i, 2 , 3 , triggering 
the Landau pole and mpi = 1.22 x 10^® GeV is the Planck mass. The states marked with | 
and -k are constrained by EW precision tests and BBN, respectively, to lie possibly beyond the 
LHG reach. 
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Spin 

X 

Qlp 

^ decay 

dim(Gdecay) 

A^rrjGeV] 

0 

(3,1,5/3) 

5/3 


5 

> "ipi {gi) 

0 

(3, 2, 5/6) 

1/3,4/3 

X^Hqq, x^Hu^e^, 
xH^u'^d^, xHu^u'^i 

5 

> mpi {gi) 

0 

(3,2,11/6) 

4/3,7/3 

xH^u^u'^, x^Hd^e^ 

5 

5.5 X 10^^ (gi) 

0 

(3, 3, 2/3) 

-1/3,2/3,5/3 

X^H^qe^-, xHW^^i 
xHU% D^^xq^a/ 

5 

> mpi (gi) 

0 

(3, 3, 5/3) 

2/3,5/3,8/3 

X^Hqe^, 

5 

3.2 X 10^" (g^) 

0 

(3,4,1/6) 

-4/3,-1/3, 
2/3,5/3 

xH^qq, X^Hqi 

5 

> "ipi (g 2 ) 

0 

(3,4, 5/6) 

-2/3,1/3, 

4/3,7/3 

X^Hqq, xH'^q^ 

5 

> mpi (g 2 ) 

0 

( 6 , 2 , 1 / 6 ) 

-1/3,2/3 

xH^qq, X^Hu^d’^, 
X^Wd^d^, Df^X^q^a^d^ 

5 

> mpi (gi) 

0 

( 6 , 2, 5/6) 

1/3,4/3 

X^Hqq, xHu^u^, 
xH^u'^d’^^ D^Xq'^^ij.'^^ 

5 

> mpi (gi) 

0 

( 6 , 2, 7/6) 

2/3,5/3 

X^Hd^d'^ 

5 

> mpi (gi) 

0 

( 8 , 1 , 0 ) 

0 

xHqu^, xH^d^, 
di^xD‘'G^^, D^^xq^^^,q, 

xG^'^G^u, xG^^Bfj,^, 
XXXH^H 

5 

> mpi (gi) 

0 

( 8 , 1 , 1 ) 

1 

xHhu^, X^Hqd<^, 
D^^X^u^^a.d^, XXX^H^H^ 

5 

> mpi (gi) 

0 

(8,3,0) 

0,1 

xHqu'^, xH^qd\ 
xG^^’^w^u, D^xq^^^lq. 

XXXH^H 

5 

> "Tpi ( 5 - 1 ) 

0 

(8,3,1) 

0 , 1,2 

xH^qu^, X^Hqd^, XXX^H^H^ 

5 

1.0 X 10^^ (g^) 

1/2 

1/2 

1/2 

( 6 ,1,1/3) 
( 6 ,1,2/3) 
( 8 , 1 , 1 ) 

1/3 

2/3 

1 

Xcr^^'u^G^^ 

5 

5 

5 

> mpi (gi) 

> mpi (gi) 
4.0 X 10^6 (^ 1 ) 


Table 2: Same as in Table 1 but for colored states. 
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Interestingly, it turns out that such SM extensions generically possess extended accidental 
symmetries which ensure the stability of the lightest particles (LPs) in the new multiplets at the 
renormalizable level. If these are charge- and color-neutral, they can form viable dark matter 
candidates, a possibility, which has been throughly investigated in the literature [4-7]. On the 
other hand, scenarios where the lightest component of the new multiplet is charged and/or 
colored are in general constrained by cosmological observations as well as by searches for exotic 
forms of matter on Earth and in the Universe. Taking also these constraints into account, 
the hnal list of viable uncolored and colored weak representations are given in Tables 1 and 2, 
respectively, which summarize the main results of our investigation. 

The details of the above sketched construction and analysis are contained in the rest of the 
paper which is organized as follows. In Sect. 2 we exploit accidental symmetries beyond the SM 
to construct SM extensions with new degrees of freedom at the weak scale, which are completely 
transparent to indirect low-energy probes. The way the set of all possible extra states is made 
hnite is discussed in this section as well. In Sect. 3 we estimate the new particles’ lifetimes. 
In turn in Sect. 4 we consider bounds on possibly long lived states coming from early Universe 
cosmology. In particular, the effects on big bang nucleosynthesis (BBN) turn out to be the 
most important ones. Sect. 5 explores the collider phenomenology of the viable weak-scale SM 
extensions and estimates the current lower bounds on the new particles’ masses coming from 
existing LHC searches. We conclude in Sect. 6 while a more detailed technical discussion of the 
renormalization group (RG) evolution of the gauge couplings and the SU(2 )l decomposition of 
the effective operators are relegated to the Appendices. 


2 Accidentally safe extensions of the SM 

Our starting point is the classihcation of SM extensions which automatically preserve the acci¬ 
dental and approximate symmetry structure of the SM without imposing additional protective 
mechanisms (only SM gauge and Lorentz symmetries are required). For simplicity, we will 
limit our discussion to the case where a single extra representation y is added to the SM held 
content. While simultaneously adding more than one representation from our set is in principle 
possible, two additional restrictions need to be considered in that case: (i) adding more matter 
representations will in general lower the scale of the EFT validity (cf. Sect. 2.4), (ii) additional 
SM gauge invariants may be constructed, potentially breaking Qp and/or the new accidental 
symmetry associated with y stability at the renormalizable level. 

We start by listing all the d < 3 operators made of SM helds. If y is a fermion, we require 
that the new state does not couple to SM fermions at the renormalizable level. In this way, 
Qp is automatically preserved and an extra accidental symmetry guarantees the stability of the 
new particle at the renormalizable level. On the other hand, the case of extra scalars is more 
involved since they can always couple to the Higgs held at the renormalizable level without 
breaking Qp and their stability depends on the allowed interactions with the Higgs held. 

A brief comment regarding larger Lorentz group representations is in order at this point. 
The presence of extra Lorentz vectors (i.e. spin 1 bosons) requires either the extension of the 
SM gauge group or new strong dynamics. In the former case, accidental preservation oi Qp 
requires the extended gauge symmetry to be a direct product of the SM gauge group and 
possible new factors under which the SM fermions need to transform trivially. Such setups 
have been thoroughly studied in the literature (c.f. [8] for a recent review) and we have nothing 
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Spin 

SM field 

SU(3), 

SU(2)i 

U(l)v- 

0 

H 

1 

2 

+1/2 

1/2 

q 

3 

2 

+1/6 

1/2 


3 

1 

-2/3 

1/2 


3 

1 

+1/3 

1/2 

i 

1 

2 

-1/2 

1/2 


1 

1 

+1 


Table 3: SM field content and quantum numbers. 



C^SM 

SU(3)e 

SU(2)i 

U(l)r 


qH{W) 

3 

1©3 

+2/3(-l/3) 



3 

2 

-l/6(-7/6) 



3 

2 

+5/6(-l/6) 


m{W) 

1 

1©3 

0(-l) 


e^H{W) 

1 

2 

+3/2(+l/2) 


qq 

3©6 

1©3 

+1/3 


qu^ 

1©8 

2 

-1/2 


qd^ 

1©8 

2 

+1/2 


qi 

3 

1©3 

-1/3 


qe^ 

3 

2 

+7/6 


u‘^u^ 

3©6 

1 

-4/3 


u^d'^ 

3©6 

1 

-1/3 

t/sMt/sM 

u^i 

3 

2 

-7/6 



3 

1 

+1/3 


d^d'^ 

3©6 

1 

+2/3 


d^i 

3 

2 

-1/6 


d^e^ 

3 

1 

+4/3 


a 

1 

1©3 

-1 


ie^ 

1 

2 

+1/2 


e'^e^ 

1 

1 

+2 


HH 

1 

3 

+1 

H 

HW 

1 

1©3 

0 

combinations 

HHH 

1 

4 

+3/2 


HHH^ 

1 

2©4 

+1/2 


Table 4: List of all possible d < 3 operators made of SM fields. Operators of the type 

are not displayed since they couple to Lorentz vectors, which are not considered in our analysis. 
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to add. On the other hand, new vectors due to some strong dynamics at the TeV scale are 
incompatible with a large mass gap Aefr S> TeV. The same argument applies to composite 
particles of higher spins. Finally, extra fundamental particles with spins 3/2 and 2 can appear 
in theories of extended and gauged space-time symmetry (c.f. [9, 10]), but such constructions 
necessarily go beyond our EFT framework. We will hence limit our discussion to the inclusion 
of either spin 0 or 1/2 extra representations. 

In the following, we adopt a two-component notation where all the fermion helds are Weyl 
spinors belonging to the same irreducible representation of the Lorentz group. The accidental 
matter multiplets are collectively denoted by x- We use the subscripts S and F to denote 
the bosonic (spin 0) and fermionic (spin 1/2) SM gauge representations, respectively, where 
appropriate to avoid ambiguity. The SM fermions are collectively denoted by V’sm and their 
quantum numbers are hxed according to Table 3. The list of all possible d < 3 operators made 
of SM helds is provided in Table 4. 

2.1 New fermions 

If a fermionic y transforms under a complex or pseudoreal representation of the gauge group 
(so that a Majorana mass term is forbidden), we introduce another held x'^ with conjugate 
quantum numbers. In this way, the new state is vector-like and a mass term can always be 
added. 

According to our previous discussion, we want to forbid the interactions xi^SM, xi^suH 
and By inspecting Table 4 we conclude that x cannot have the following quantum 

numbers: 

X ^ ^SM, (1,1, 0), (1, 3, 0), (1, 3,1), (1, 2, 3/2), (3, 2,5/6), (3, 2, 7/6), (3, 3,1/3), (3, 3, 2/3). (2) 

If X transforms under a real representations of the SM group, then we can also add a Majorana 
mass term and the most general Langragian reads (see e.g. [11] for two-component notation) 

c = CsM + ix^a^D^x + ^^(x^W + h-c.), (3) 

which is invariant under a Z 2 transformation x —x- On the other hand, if x transforms 
under a complex or pseudoreal representations of the SM group, we introduce an extra Weyl 
fermion x^ with conjugate gauge quantum numbers with respect to X) so that a Dirac mass 
term is allowed, and get 

^ = ^SM + ix^a^^D^x + + Mix^ex"" + h.c.), (4) 

which is invariant under a U(l) transformation x e*^x ^md x'^ both cases an 

accidental symmetry implies stability of the new particles at the renormalizable level and also 
requires that they are pair produced in high-energy particle colliders. 


^Notice that terms of the form xxH or xxH^ s-re forbidden by SU(2)j^ invariance. 



2.2 New scalars 


For scalar x, in order to preserve Qf we have to avoid all couplings of the form xV’smV'sm- By 
inspecting Table 4 we conclude that y cannot have the following quantum numbers: 




( 1 , 1 , 

(3,3, 


1),(1,3, 
1/3), (6, 


1), (1,1,2^), (1, 2,1/2), (3,1,1/3), (3,1, 2/3), (3,1,4/3), (3, 2,1/6), (3, 2, 7/6), 
1,1/3), (6,1, 2/3), (6,1,4/3), (6, 3,1/3), (8, 2,1/2). (5) 


Spin 

X 

^decay 

dim(C>decay) 

Stability 

0 

(1,1,0) 

XHW 

3 

X 

0 

(1,3,0) 

xhW 

3 

X 

0 

(1.4,1/2) 

xHWW 

4 

X 

0 

(1,4,3/2) 

XH^H^H^ 

4 

X 

0 

{R,2k,l/2) 

XXH^H^ 

4 

Z2 

0 

(i?,n,0) 

XXHH^ 

4 

Z2 

0 

(C, n, Y) 

xx^hW 

4 

U(l) 

0 

(C',2fc,l/6) 

XXXH^ 

4 

^3 

0 

{R,2k,l/2) 

XXX^H^ 

4 

X 


Table 5: Extra scalar representations which can couple to the Higgs at the renormalizable level 
without breaking Qp- {C,n,Y) denote generic quantum numbers under the SM gauge group 
which are not already contained in the list of Eq. (5). R stands for a real SU(3)c representation 
(i.e. R = 1,8, 27,...) and 2k for an even 811(2)^ representation. In the last column, we provide 
(when appropriate) the symmetry responsible for the stability of y. The cases denoted by a 
“x” lead instead to the decay of y at the renormalizable level. 

Analogously to the case of extra fermions in Sect. 2.1, gauge interactions alone cannot lead 
to the decay of x at the renormalizable level, since the kinetic terms again exhibit a Z 2 or a 
U(l) invariance for the case of an extra real or complex scalar, respectively. The decay of the 
new particle is however possible (depending on the quantum numbers of x) due to the presence 
of extra renormalizable interactions between x and H, which are listed in Table 5. 


2.2.1 Scalar potential, CP and custodial symmetry 

In the presence of any new scalar multiplet x the scalar potential can be written as (see e.g. [12]) 

V{H,x) = VsM + <1 K 1^1" + “ Y + ^(■X'T°x){H'Ti,H)) 

+ [x(x'C-,T‘x‘){H'^C„TiH) + hx.] + ... , (6) 

where rj is equal to 1(1/2) for a complex (real) representation, and Cr denote respectively the 
SU(2)i generators and conjugation matrices in the representation R (so, for instance, = cr“/2 
and Ch = where a°‘ for a = 1, 2,3 are the Pauli matrices). We take (H)'^ = (0, with 

V = 246 GeV. The ellipses in Eq. (6) stand for extra terms, like e.g. {x^Tl^xYi which do not 
sizeably affect the mass splitting of x (see below). In addition, x in specihc weak representations 
might allow for additional renormalizable operators listed in Table 5. 

The hrst accidental symmetry of the scalar potential that we wish to discuss is CP. Generic 
sources of CP violation are severely constrained by the measurement of electric dipole moments 
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(EDMs) [13]. Among the accidental scalar matter extensions of Tables 1-2, it turns out that 
only (1,4, 1 / 2)5 explicitly violates CP. This can be seen by noticing that for such a multiplet one 
can construct three non-hermitian invariants in the scalar potential (cf. the third, hfth and last 
row in Table 5) and that only one out of the three phases associated with the corresponding 
complex couplings can be rotated away by a re-phasing of x and H. In this case the most 
signihcant experimental constraint comes from the searches for an electron EDM (dg), dehned 
through the effective operator C 3 —i{de/2)e{a ■ F) 75 e. The (1,4, 1 / 2)5 contributes at two 
loops through the diagram in Fig. 1, which corresponds to the diagram in Fig. 12 of Ref. [14] 
after replacing Hi ^ H and H 2 -3 x- 


W+' 


/ \ 

/ \ 



X+^ 

^X+ 

/ 

\ 

/ 

^ Xo 



Xo 

Ro A 




\ 

/' 



Figure 1: Scalar loop contribution to the electron EDM. 


Taking into account the extra (y) /u < 1% suppression due to EW precision constraints 
(see below), assuming 0{1) scalar couplings and mixing angles, and ~ v for the sake of 
a very conservative estimate (see also [14, 15]), we obtain |de| < 7 x 10“^®ecm. This has 
to be compared with the recent experimental bound from the ACME collaboration [16] of 
l^expj ^ ^ 10“^®ecm at 90% C.L.. While not constraining at the moment, interestingly, 

future experimental improvements on the electron EDM might start to probe CP violation in 
generic weak-scale scalar extensions of the SM involving the (1,4, 1 / 2)5 multiplet. 

Another accidental symmetry of the SM scalar potential is the so-called custodial symmetry. 
In the 5 ^' —)■ 0 limit the massive gauge bosons transform as a triplet of an unbroken global 
811 ( 2 ) 0 , which is also responsible for the tree-level relation ptree = m^/m|cos^ 6 *w = 1. New 
sources of SU(2)o breaking which cannot be accounted in the SM are described by the po = 
p/psM parameter [17]. Experimentally, = 1.0004/Q;oQg4 [17], which is compatible with the 
SM prediction po = 1- Thus the experimental value of po can be used to constrain new sources 
of SU(2)o breaking due to the extra scalar x- 

If X gets a vacuum expectation value (VEV), there is a tree-level contribution [18] 


pr 


l = {p[j{j + l)-Y^]-2Y^] 


{xY 


o 


{xY 


(7) 


where j is the total weak-isospin quantum number of x cind Y its hypercharge in the Q = 
+ Y normalization. Apart for the safe representations yielding po = 1 for any value of (x): 
(1,1,0), (1,2,1/2), (1,7,2) [19], (1,26,15/2), etc., the 2cT-level saturated bound is at the level 

of (x) h < !%• 
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In general, whether a scalar field can develop a VEV depends on the choice of the parameters 
in the scalar potential. However, “tadpole” couplings of x to some hf’s always imply an induced 
VEV for X- From Table 5 we see that this is indeed the case for the states: (1,1,0), (1,3,0), 
(1,4,1/2), (1,4, 3/2). While the VEV of the former does not contribute to po, the remaining 
ones can be in principle dangerous. By looking at the generic shape of the potential and its 
stationary equations, we estimate on dimensional grounds (for 0{1) couplings and barring hne- 
tunings), (x) ~ v^/m^ (triplet case) and (x) ~ v^/rn^ (quadruplet cases). Hence, (x) /v ^ 1% 
corresponds to > 100 v ^ 20 TeV (triplet case) and > 10 n ~ 2 TeV (quadruplet cases), 
which limits the visibility of these states at the LHC, unless a moderate hne-tuning is allowed 
in the scalar potential. 

Custodial symmetry also helps us to understand the properties of the theory beyond the 
tree level. Indeed, a tree-level splitting within the components of x originating from the scalar 
potential in Eq. (6) gives a radiative contribution to po- In the following we assume (x) <C f to 
suppress the tree-level contribution to po- Consequently (x) itself cannot sizably contribute to 
the mass splitting. Notice, also, that among the scalar states selected in Table 1, the coupling 
7 is relevant only for (1,4,1/2). However, since this state decays through a renormalizable 
operator, the details of its mass spectrum are not of particular interest.^ We are hence left with 
the contribution of (3 to the mass splitting, which yields 

mj = — 5 ^/, ( 8 ) 

where —j<I< j denotes the eigenvalue of the ( 2 j -|- l)-dimensional representation x and 
we defined the parameter -|- ^av'^ and Using the general formula for the 

one-loop correction in [20] and expanding the loop function for 6 < M we find 


1—loop 

Po 


- 1 


pNcO-era 2(5 Y'|iVO'i1'\iCO 



(9) 


where Nc is the dimensionality of x under the color factor. Neglecting the higher-order 5/M 
terms, we hnally obtain 


M > 72.5 GeV 


0.001 

exp ^ 

Po - ■ 


1/2 


^VvNcVjU + l)(2j + 1), 


( 10 ) 


which is valid for M > 5 ^ y/]3 123 GeV. For 0{1) values of the coupling f3 the typical 
bounds on M range in the few hundred GeV region, depending on the dimensionality of the 
representation. We hence conclude that the mass bounds coming from loop-level contributions 
to Po are less general (they depend on the value of f3) and not particularly constraining when 
compared to existing direct searches limits (see Sect. 5). This is, however, not necessarily true 
for higher dimensional representations. 

Alternatively, the information from po can be used to give an upper bound on the mass 
splitting Am = mj^i — mj ~ As an example, let us mention that for the case (1, 5, 2)5 

we get Am < 20 GeV. This information is exploited in Sect. 5.3 when inferring collider bounds 
on the neutral state of such a multiplet by looking at the charged component production and 
decays. 

^ The coupling 7 induces mixing between the conjugate components of x with the same \Q\ A 0 and, for the 
Q = 0 component, it splits its real and imaginary part. The contribution of 7 to has been considered 

for instance in [ 12 ]. 
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2.2.2 Bounds on Higgs portal coupling 

The Higgs boson can couple to the new scalars via the portal coupling a of Eq. ( 6 ). This 
leads to two kinds of effects: (1) If some components of x lie below half of the Higgs mass, 
they can contribute to the Higgs total decay width. Taking into account other existing collider 
constraints (see Table 8 ) this is only possible for the neutral component Xo- Iii particular, it 
contributes to the Higgs invisible decay branching fraction. The partial decay width of the 
Higgs boson via the a coupling into a pair of xo states (for /3 = 0) is found to be 


IGtt ruH y 


( 11 ) 


In addition, (2) all charged components of x will contribute at 1-loop level to the if —)■ 77 (and 
H —)■ 7 Z) decays, while colored x will affect Higgs boson production through gluon fusion (GF) 
and also its decays to two gluons H ^ gg. Using the results of [21-23] we hud 


_ T^^ _ \Ai{xw) + {4:/?,)Ai/2{xt) + gad{R^) 

'‘"“7“ |A(i„,) + (4/3)A/27.)P 

^ ^ ^ ^ \{M‘^)Ai/2{xt) + vaC{R^) T.iiv/mxi)^Ao{x^,)\'^ 

'‘"“7“ a*' i(i/2)A/2(i.)r^ 


where Xi = the sums over all x weak multiplet components Xi, d{R^ is the 

dimension of the color representation of x and C{R^) is the corresponding index (G(3) = 1/2, 
(7(6) = 5/2 and (7(8) = 3). The relevant loop functions Ai{xw) — —8.32, Ai/ 2 {xt) — 1.38 and 
Ao{x) with limits Aq{x —)■ 0 ) = 1/3, Ao{x —)■ 00 ) = —1/x + 0(x~^) can be found e.g. in [23]. 
The total decay width of the Higgs can thus be written as 


Fff = [1 + - 1 ) + - 1 )] + Ti,, , (13) 

where Tf/^ = 4.07 MeV, = 2.28 x 10“^ and HR™ = 8.57 x 10“^ [24]. The invisible 

branching ratio is then finally given by BRinv = Tinv/T/f. 

To analyze the resulting correlated effects in Higgs boson production and decays as mea¬ 
sured at the LHC, we follow closely the procedure described in [26] using also the same set of 
experimental results [25, 27-36]. In particular, we hud that in cases where the Higgs boson can 
decay to Xoi the constraints on a are completely dominated by the bounds on the extra invisible 
decay rate. In Fig. 2 the invisible branching ratio as a function of the new scalar mass is 
shown for different values of the portal coupling a. We used 7 = 1/2 in Eq. (11), assuming 
a real scalar. The red dashed line shows the CMS limit of Ref. [25].^ It can be inferred from 
the plot that for a portal coupling |q;| = (7(1), the new scalar states are excluded up to the 
kinematic limit for this decay. However for values of |a| < (7(0.01) currently no limit on 
can be given anymore. 

Even if Xi are heavy (m^. > mn/^)^ their contributions to and figg still lead to constraints 
on a from the measurements of the Higgs signal strengths at the LHC. In particular, the most 

^Indirect bound on BRinv coming from the global fit to all Higgs boson signal strenghts yields a slightly 
stronger bound of BRinv ^ 0.2. 
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Figure 2: Invisible branching fraction of the Higgs boson as a function of the new scalars mass 
(real scalar) for different values of the portal coupling a. The red dashed line shows the CMS 
exclusion limit of Ref. [25] . 


sensitive channels involve GF produced Higgs bosons decays to photons and W bosons, these 
being the two most precisely measured. Denoting the relevant signal strenghts as 


/i 


GF _ 
77 ~ 


(Jgf hr, 




SM 

GF 


HR: 


77 
SM ’ 
77 


,,GF _ '^GF BRww 
hww = ZMl d-dSM 


0 -; 


GF 


HR: 


ww 


(14) 


the global £t of Higgs boson LHC data allowing for arbitrary contributions to and figg but 
keeping Fjnv = 0 yields the 68% and 95% CL exclusion bounds shown in Fig. 3 . We observe 
that up to 50% modihcations in both observables are still allowed by the current data. These 
should be compared with a induced modihcations shown in Fig. 4 (assuming degenerate Xi)- 
In particular, color-neutral y predominantly affect as shown in the left panel. On the 
other hand, colored states can affect GF production and are thus constrained also from 
as illustated in the right panel. The deviations in are shown for single complex scalar 

in the given color representation. Finally, colored scalar effects in are also shown in the 
middle panel. 

Asymptotically, y effects in both observables decouple as a/rri^. The shaded bands in Fig. 4 
illustrate the amount of deviations from this limit as they correspond to a scan of |q!| G [0.1,1]. 
We observe that for > 500 GeV even |a| < 0{1) can be consistent with current Higgs data. 
Gonversely > 100 GeV are perfectly allowed for small enough |a| < 0(0.1) Higgs portal 
couplings. 


2.3 Mass spectrum 

The phenomenology of the new EW states is dictated by the mass spectrum. Typically, on top 
of a common mass term there is a radiative splitting within the 811(2)^ multiplet and, for 
scalars only, a tree-level splitting due to the presence of non-trivial SU( 2 ) 2 , invariants in the 
scalar potential. In the 3> v limit the radiative contribution takes the form [4, 37] 

Amrad = mq+i - tuq Ri 166 MeV (l + 2Q+ ] , (15) 

V cos Ow J 
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Figure 3: Exclusion bounds on the and LHC Higgs signal strengths allowing for 

arbitrary contributions to and figg but keeping Finv = 0. 




m;j,/V I a I [GeV] \ a \ [GeV] m^/V | a \ [GeV] 


Figure 4: Accidental scalar induced modihcations to and LHC Higgs signal strengths 
as a function of m^/ ^/\a\. The un-colored and colored scalar effects in are shown in the 
left and middle panel, respectively. Single (complex) colored scalar effects in are shown 

in the right panel. The shaded bands correspond to a scan |a| G [0.1,1]. 


which holds both for fermions and scalars. Notice that if H = 0 the LP in the multiplet is 
always the one with the smallest |Q|. This is not necessarily true when Y ^ 0.^ 

Similarly, the tree-level splitting in Eq. (8) can be expanded in the S> v limit, thus 
obtaining [4] 

Bv‘^ fl TeV\ 

Amtree = - mj ^ ^ /3 X 7 .6 GeV ( - ) . (16) 

Notice that, while for fermions the mass spectrum is unambiguously hxed, for scalars it depends 
on the values of (3 and m-^. Focussing on the < TeV region (relevant for LHC), if B = 0{1) 
then the LP is always the one with the highest/lowest /, depending on the sign of B- However, 
for B £ [10”^) 1] fhe tree-level splitting can be comparable with the radiative one. In such cases 
it is possible to show (see below) that: 

1. Any particle in the multiplet can be the LP for large domains of the model parameters, 
i.e. without any hne-tuning. 

®E.g. the LP of the fermion multiplet (1,4, l/2)p has Q = —1. 
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2. If the LP has charge Qlp, the next-to-LP has always charge Qlp i 1. 

This latter fact turns out to be phenomenologically relevant, e.g. when setting bounds on the 
neutral LP by looking at the decay of the next-to-LP. 

For completeness, we provide here a proof of the two statements above: by combining 
Eq. (15) and Eq. (16) one arrives at the expression mj = rri-j + a + bl + cP, where a and b can 
have any sign (since they depend on Amtree), and c > 0. The minimum of mj is obtained for 
Lmin = —Hence, by an appropriate choice of the ratio b/c, the smallest mj can be anywhere 
in the range I G The fact that the next-to-LP has Qlp =t 1 simply follows from the 

convexity of mj as a function of I. A similar argument holds as well in the ~ v regime, 
for which the full formula of the radiative splitting (see e.g. Eq. ( 6 ) in [4]) must be taken into 
account. 

2.4 Validity of the EFT 

Our working hypothesis is that the SM-l-y renormalizable theory is a low-energy effective de¬ 
scription valid up to a cut-off scale Aeg. In the spirit of a generic EFT with 0{1) couplings 
and without any extra state beyond y introduced at low energy, Agff ~ 10^® GeV is essentially 
hxed by neutrino masses through the d = 5 Weinberg operator. Moreover, such a cut-off scale 
can automatically account for null results of all flavor, CP and B violating processes constrain¬ 
ing d = 6 operators made of SM helds. In particular, when the lowest-dimensional sources of 
breaking of the extra U(l) or Z 2 symmetry associated with the kinetic term of y are the d = 5 
operators involving y and SM fields, any d = 6 operator involving only SM helds, generated by 
integrating out y, will have two insertions of such d = 5 operators and hence at least a l/Ag^ 
suppression. The situation changes only slightly if the extra U(l) or Z 2 is broken at the renor- 
malizable level in the scalar potential, as in cases listed in Table 5. Namely, the only additional 
effect arises for y ~ (1,4, 1 / 2 ) 5 , (1,4, 8 / 2)5 where integrating out y induces a AL = 2 operator 
of the form UHHHW. Being suppressed by l/Aegm^, it necessarily represents a sub leading 
contribution to neutrino masses. 

The inhnite set of states preserving Qp the renormalizable level (see Eq. (2) and Eq. (5)) 
can be reduced by requiring that the EFT remains weakly coupled up to Agff ~ 10^^ GeV. 
The presence of extra matter multiplets drives the gauge couplings of the SM towards the non- 
perturbative regime.® Eventually, this might result in the presence of a Landau pole below the 
cut-off scale of the EFT. If the Landau pole is associated with a generic new dynamics, the 
accidental symmetries of the SM could be violated at that scale. Hence, for the self-consistency 
of the EFT approach, we require the absence of Landau poles below Agfj 10^® GeV, which 
translates into an upper bound on the dimensionality of the extra representations. 

In light of stringent bounds on the inter-multiplet mass splittings (see Sect. 2.3) we can 
safely integrate in all multiplet components at a single scale, which we choose to be the Z mass 
in our numerical analysis. We note however, that for not much larger than the TeV scale the 
resulting Landau pole estimates scale linearly with y masses. The analysis of the perturbativity 
bounds is detailed in Appendix A and the results are summarized in Tables 9-10. They provide 
a useful reference for the estimate of the Landau poles at two loops for the cases where the SM 

®We do not address here the question of the RG running of the scalar potential parameters, since it is a 
model dependent issue which also involves the analysis of the vacuum stability. 
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is extended with an extra multiplet charged under SU(3)c and/or 811(2)/,, and in particular for 
all the states considered in this work which can have a non-zero hypercharge as well. 

A crucial ingredient in order to make our list of extra states hnite however, is given by 
cosmology. In fact, the only reason why we can disregard multiplets with an arbitrary hyper¬ 
charge, e.g. F = TT, is because these states feature an absolutely stable charged LP that cannot 
decay into SM particles because of electric charge conservation. The possibility of having an 
inhnitesimal hypercharge is instead briefly discussed in Sect. 4. 

A comment on the role of higher-order corrections in the RG equations is in order here. 
The determination of the Landau pole is often carried out at the one-loop level (see e.g. [4]). 
However, for the non-abelian gauge factors there is an accidental cancellation in the one-loop 
beta function between matter and gauge contributions (cf. Eq. (36)), so that two-loop effects 
may become important. Interestingly, among the cases that we found to be drastically affected 
by two-loop corrections there are the two minimal DM candidates: a real (1, 7, 0) scalar and a 
Weyl (1,5,0) fermion [4].^ Following the results of [5] for the calculation of the relic density, 
we integrate in the scalar septuple! at = 25 TeV and the fermionic quintuplet at = 10 
TeV. Hence we hnd, respectively 

A tS ,, = 1.9 X 10^‘ GeV ^ Afi = 8,9 X 10"" GeV ((1,7,0)s case) , (17) 

ASL = 9,0 X 10"" GeV ^ a S,, = 4.0 x 10"‘ GeV ((1,5,0)js case) . (18) 

If we associate the Landau pole with the cut-off of a generic EFT, this also sets the scale of the 

effective operator leading to the decay of the minimal dark matter candidate. Note however, 

that even for a cut-off of the order of the Planck mass, the framework of minimal DM is not 
endangered by d > 6 operators, since the lifetime of DM is still comfortably larger than the age 
of the Universe (and satishes the indirect bounds on decaying DM). For a discussion of d = 5 
induced (1,7,0)^ decays see Sect. 3.3.2. 

In the selection of our states, the two-loop criterium proved to be important for several 
states. For instance, in the case of the real (27,1, 0)5 scalar multiplet we hnd that at one loop 
^Landau ~ ^ GeV, whereas at the two-loop level ALa°d^^ = 1.3 x 10^ GeV, so that we 

can exclude this state from our list of accidental matter candidates. 

What about three-loop corrections then? As long as there are no accidental cancellations in 
the two-loop beta function (as it can be explicitly verihed), they are not expected to drastically 
change the situation.® It is then enough to rely on a two-loop estimate of the Landau pole in 
order to set an upper bound on the dimensionality of the extra representation. 

3 Lifetimes 

The new extra states will eventually decay due to operators present in the EFT. There are 
essentially three classes of decays which we are going to consider in this section: i) Inter- 
multiplet weak transitions where the heavier components within the SU(2)/, multiplet decay 

^Another situation where the two-loop RG analysis of the gauge couplings could change the qualitative UV 
behaviour of the theory is given by the Pati-Salam model presented in [38], where low-scale extensions of the 
SM providing total asymptotic freedom are investigated. 

®For instance, in the SM case where no strong cancellations are at play we find: = 1.9 x 10^^ GeV, 

= 5.2 X 1040 GeV and A?±°p , = 8.7 x lO^o GeV. 
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via cascades involving the emission of (virtual) W gauge bosons into the LP, iij Decays through 
renormalizable interactions (only for a specihc class of new scalar states) and iiij Decays through 
non-renormalizable d > 5 operators. We analyze each class of decays in turn below. 


3.1 Inter-multiplet weak transitions 

Heavier components within the 811(2)^ multiplet can decay via EW transitions into lighter 
ones, with rates suppressed by a small phase space factor. Denoting the component of a total 
j-isospin representation with T^-eigenvalue I as xh we have the decay width 

(generalizing the expression in Ref. [37]) 


r(xP. ^ x‘, ir+) = 




TT 


1 - 


m 


Am2 7.5 X 10 


-12 


Am \ ^ 
500 MeV ) 


(19) 


where T+ = xj]{] + 1) — /(/ + 1) and the approximation in the r.h.s. of Eq. (19) is valid for 
Am m 7 r+. 

Formula (19) is a reasonable approximation of the total width in the range mjr+ < Am < 
1 GeV. For mass splittings close to the kinematical threshold of the decay into a pion, 3-body 
decays involving leptons become important as well, while for Am > 1 GeV new hadronic chan¬ 
nels open up (e.g. involving kaons and other heavier hadrons) and the decay can be eventually 
computed at the partonic level, once quark-hadron duality sets in. 

The typical lifetime of an SU(2)/, multiplet component decaying via inter-multiplet weak 
transitions is displayed in Fig. 5 as a function of the mass splitting and for different values of 
the ladder operator T+, up to the j = 3 (septuplet) case. 

Within high energy collider experiments, the inter-multiplet decays are essentially prompt. 
On the other hand, the LP at the end of these inter-multiplet cascades is stable on the detector 
scale, barring few exception which are discussed in the next subsection. 


3.2 Decays through renormalizable interactions 

There exists the possibility that the new extra scalars retain renormalizable interactions with 
the SM Higgs which can induce their decay. These states are classihed in Table 5 and correspond 
to the cases (labelled with the symbol “x”) where no accidental symmetry (e.g. U(l), Z 2 or 
Z 3 ) forbids x to decay. Let us comment in turn on the various possibilities. 

The case of the gauge singlet (1,1, 0)5 has been extensively studied in the literature (see 
e.g. [39]) and we do not have much to add here. In the cases (1,3, 0 ) 5 , (1,4,1/2)^, and 
( 1 , 4 , 372 ) 5 , X can couple linearly to Higgs operators. However, such “tadpole” couplings also 
induce non-zero VEVs for x, which are severely constrained by EW precision observables. As 
already pointed out in Sect. 2.2.1, unless a moderate hne-tuning is allowed in the scalar poten¬ 
tial, the bounds on such dimensionally estimated VEVs push the masses of these states beyond 
the kinematic reach of the LHG. Other multiplets which can possibly decay at the renormal¬ 
izable level (those labelled with a “x” in Table 5) are either not considered here because they 
break Qjr by coupling to SM fermions (e.g. (8,2, 1 / 2 ) 5 ) or they generate a Landau pole below 
GeV (e.g. ( 1 , 6 , 172 ) 5 ). 
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Figure 5: Lifetimes associated with inter-multiplet weak transitions as a function of the mass 
splitting Am. The grading of the curves (from black to gray) corresponds to different values 
of (from 1 to 12), as it can be found in representations up to j = 3 (septuplet). The red 
dashed line corresponds to the typical freeze-out time for colorless 1 TeV-mass particles with 
weak interactions in the early Universe (cf. Sect. 4). 

3.3 Decays through d > b effective operators 

Let us consider now the case where the decay of the new state y is due to effective operators. 
Given an effective operator Odecay, we always absorb the Wilson coefficient in the dehnition of 
the effective cut-off scale Agff, e.g. 


C 3 


1 


A, 


-a 


eff 


decay 


h.c.. 


( 20 ) 


The differential decay rate of an unstable particle x into nf hnal states reads 

^ ( n ^ (27r)V^) (p^ 


2 m, 


/ 


/ 


( 21 ) 


/ 


By assuming a constant matrix element and massless hnal states, the phase space factor can 
be integrated in the rest frame of the decaying particle, yielding 
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d^Pf 1 


27r3 2E 


f 


{271^5^^^ [px-J2Pf 
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2(47r)2»i/ ^ — 2)\ 


( 22 ) 


So, for example, the phase space factors up to rij = 4 are: PS 2 = PS 3 = ^^^ 3 , PS 4 = 24 ^^- 
In the case of a dimension d effective operator, the amplitude squared for Uf particles in 
the hnal state can be estimated by naive dimensional analysis (NDA) as 


|A4(m^ -3 {p/})Inda 


/ V \2nr. 

'\/2'^ 2d—2nf—2nc—6 
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.2(d-4) X 


(23) 
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where we also included the possibility of ric condensations of the Higgs boson. Hence, by putting 
Eqs. (21)-(23) together, we get the following expression for the total width for 3> v 


Tnda — 


m 


2d-2nc-7 (jv_\2nc 
X 


4(47r)2”-/ 3 — 2)! 


(24) 


Unless differently specified, we compute the lifetimes of the states decaying through the non- 
renormalizable operators in Tables 1-2 using Eq. (24). Whenever multiple operators can be 
responsible for the decay of y, we sum over the several widths assuming the operators contribute 
with the same Wilson coefficient. What is missing in Eq. (24) with respect to the full decay 
width are the relevant SU(2)i Clebsch-Gordan coefficients, symmetry, color and flavor factors, 
the kinematical dependence of the matrix element, the masses of the hnal states and hnally 
mixing effects induced when scalar y obtain VEVs. In the region 3> v all of these are 
expected to give 0{1) corrections. When more accuracy is required, for example when setting 
BBN bounds, we take all these factors into account, computing the relevant decay widths 
explicitly. 


3.3.1 Cascade decays 

Whenever a Higgs doublet is contained in a SM-invariant operator, it can happen that not all 
the SU(2)i components of the multiplet y can directly decay through the effective operator. 
This is easily understood by going to the unitary gauge, where some of the SU(2)i contractions 
end up into the goldstone directions of the Higgs doublet. See Appendix B for a description 
of the SU(2)j;, decompositions of the relevant operators. Depending on the mass spectrum, 
the cases where the LP cannot directly decay through the effective operator are displayed in 
Table 6. 


Spin 

X 

Qlp 

C^decay 

0 

(1,2, 5/2) 

3 

X^He^e^ 

0 

(1,5,1) 

-1,1,2,3 

x^hhhW 

0 

(1,5,2) 

1,2, 3,4 

X^HHHH 

0 

(3,2,11/6) 

7/3 


0 

(3, 3, 5/3) 

8/3 

X^Hqe'^ + xH^W^^ 

0 

(3,4,1/6) 

5/3 

xH^q + X^Hqi 

0 

(3,4, 5/6) 

7/3 

X^Hqq + xHU^ 

0 

(6, 2, 7/6) 

5/3 

X^Hd^d^ 

0 

(8,3,1) 

2 

xH^qu^ + x^Hqd’^ 


Table 6: Extra states decaying through off-shell cascades. 

It is possible, however, for the LP to cascade decay via off-shell heavier components (which 
eventually decay through the effective operator) and W bosons, thus resulting in lifetimes which 
are typically larger than in the case of the direct decay. Moreover, these decay rates must be 
evaluated numerically since the NDA formula in Eq. (24) cannot be straightforwardly applied 
due to the strong momentum dependence of the matrix element. 

For the computation of the decay width of the (1, 2, 5 / 2)5 multiplet component with Q = 3, 
X+3, under the assumption that x+3 is lighter than (or degenerate with) x+2, we take into 
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Figure 6: Feynman diagrams for the cascade decay of the x+3 component of the (1,2, 5 / 2)5 
state for 


account the decays into two leptons and into two leptons together with a Higgs boson. The 
relevant Feynman diagrams are shown in Fig. 6 . The numerical phase space integration is 
performed with the help of RAMBO [40] and we neglect the effects of lepton masses. 

For the (1,5, 1)5 and (1,5, 2)5 multiplets, longer decay chains are possible for the multiple 
charged components of the multiplet. In addition, for the neutral states within these multiplets 
the effective operator induces a mixing with the Higgs boson which in turn generates direct 
couplings to SM vector bosons. These contributions to the lifetimes do not decouple for large 
and hence need to be taken into account over the whole considered mass range. For all multiplet 
components we thus consider decays with hnal states comprising of 2 — 4 SM gauge or Higgs 
bosons. The numerical results have been obtained with Madgraph 5 [41] using FeynRules [42] 
generated model hies. 

Finally, the cascade decays of the colored cases can be estimated from the one of the x +3 
component of the (1,2, 5 / 2)5 multiplet by appropriate replacements of Clebsch-Gordan coeffi¬ 
cients and by multiplying with the respective color factors. Note that an accurate evaluation 
of the decay rates for the colored cases is not necessary since their relic abundance turns out 
to be very suppressed resulting in no relevant BBN constraints. More details can be found in 
Sect. 4. 

The 811(2)+ factors needed in the evaluation of the cascade decays are exemplihed in Ap¬ 
pendix B. For all the cases the cut-off scale Aeff was set to 10^® GeV. We do not include off-shell 
effects of the W bosons in the computation of the lifetimes.® 

In Fig. 7 the lifetime of the Q = 3 component of (1,2, 5 / 2)5 is shown as a function of its 
mass, assuming that it is lighter than the Q = 2 component and hence decays via an off-shell 
X+ 2 - The blue dashed curve shows the lifetime in the zero mass splitting approximation, while 
the red solid one stands for a NDA estimated tree-level mass splitting as given in Eq. (16). In 
the plot we assume that the d = 5 operator involves only one lepton flavor. If x couples in the 
same way to all three flavors the corresponding lifetimes are reduced by a factor of three. From 
Fig. 7 it can be inferred that the presence of tree-level mass splitting only affects the lifetimes 
for low masses of x+ 3 ; of the order 0(v). For larger masses it quickly becomes completely 
irrelevant and we henceforth work in the zero mass splitting limit whenever computing cascade 
decays. 

In Fig. 8 , we show the lifetimes of all the components of the (1, 5, 1)5 (left panel) and (1, 5, 2)5 
(right panel) multiplet assuming inter-multiplet mass degeneracy. As it can be inferred from 

®In Ref. [43] it was shown that in the case of stop decays these off-shell effects can be numerically relevant 
for mass differences between the decaying particle and the decay products up to 35 GeV. 
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(1,2, 5/2)s, Q=3 



Figure 7: Lifetime as a function of the mass of the Q = 3 component of the multiplet (1, 2, 5/2), 
if it is the lightest. The red solid curve shows the lifetime for a NDA estimated tree-level mass 
splitting between x +3 and x +2 components, while the blue dashed curve represents the zero 
mass splitting limit. 


the plot, the same scaling behavior of all the components for large is found as expected in 
the SU( 2 ) 2 , limit. The lifetimes of x+2, X+3 (and x+4 in case of (1,5, 2 ) 5 ) are larger due to the 
fact that a smaller number of hnal states is available, especially at lower masses, and hence the 
decay widths are suppressed. For such long lifetimes there are potential issues with cosmology 
(see Sect. 4.2). 

3.3.2 Loop-induced decays 

In all the SM extensions considered in Tables 1-2 there is always an operator responsible for the 
decay of the new multiplet that is linear in y, except in the case of the (real) scalar multiplet 
with SM gauge quantum numbers (1,7,0). In this case the operator responsible for y decay 
is xXXH'^H}^ This can be understood by simply noticing that the SM extended with a real 
(1,7,0) scalar has an accidental Z 2 symmetry, y —)■ —y, at the renormalizable level and the 
presence of an operator trilinear in y clearly breaks such a symmetry. We note that in the 
context of minimal DM [4] this d = 5 operator and its effect on the scalar septuple! lifetime 
have been previously overlooked. The decay only proceeds at one-loop level and, depending on 
the nature of the lightest particle in the multiplet, can result in the following hnal states with 
EW gauge bosons^^ 

• yo: the possible two-body hnal states are 77 , 7 Z, ZZ and W~^W~. The relevant Feynman 
diagrams are shown in Fig. 9. By neglecting the gauge boson masses in the hnal state we 

^^Different SU(2)i contractions give rise to different independent operators. In this section we consider the 
case where two fields y are contracted in a j = 4 weak isospin multiplet. 

^^For very large y masses, final states containing Higgs bosons might be important as well. 
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Figure 8 : Left: Lifetimes of the Q = 0 (dark blue dashed), Q = I (yellow large dots), Q = —1 
(turquoise solid), Q = 2 (pink dash-dotted) and Q = 3 (violet small dots) states of the (1, 5, 1)5 
multiplet, assuming for each curve that the respective component is the lightest one. Right: 
Same as for the left-hand side but for (1,5, 2 ) 5 . The turquoise solid line corresponds to the 
Q = A lightest state. 


Xo 


/VWW 


\ / 
- - X 


\ 


XAAAAy y 



Figure 9: Feynman diagrams for the loop decay of the xo component of the (1, 7, 0)5 multiplet 
with VV = 77 , jZ, ZZ, W~^W~. Electroweak VEV insertions are denoted by a cross. 
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where Cq ~ —0.0966 is a numerical factor coming from the evaluation of the relevant 
Passarino-Veltman functions. We observe that even for an EFT cut-off at the Planck 
scale, the fast decay of the neutral component of the septuplet effectively rules out this 
particular minimal scalar DM candidate [44]. 

X+i- the two-body hnal states are IF 7 or WZ, with a decay rate given by 
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X+ 2 - there is only a two-body decay into WW, yielding 
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(27) 
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• x+ 3 - ill this case there are no two-body decay channels into gauge bosons, while it is 
possible to show that if we ignore the effect of SM fermions x+s cannot decay into an 
odd number of gauge bosonsd^ Hence, we do expect that the leading contribution to this 
decay will be given by a hnal state containing four gauge bosons. Though we did not 
explicitly compute this decay rate, we can quote (and use in the numerical analysis) an 
NDA estimate given by 


, ^ 1 

31457287r7 


1.9 X 10"^^ s“^ 


/ lO^^ GeV 

V ^eff 


/ lTeV \ 

V "lx / ’ 


(28) 


We end this section by noting that for Aefr > 10^® GeV the loop-induced lifetimes when com¬ 
bined with cosmological considerations preclude the scalar septuplet to be within the kinemat- 
ical reach of the LHC (see Sect. 4.2 for details). 


4 Cosmology 

Most of the accidentally safe SM extensions are characterized by color- and weak multiplets 
of scalars or fermions, with weak-scale masses and no renormalizable interactions beyond their 
couplings to the SM gauge bosons (and the Higgs portal operators in the case of scalars). Thus 
they will be produced and thermalized in the early Universe, eventually freezing-out once their 
thermalizing interactions become slower than the Hubble expansion rate H{T). The details 
depend somewhat on the mass hierarchy within the y multiplets but the decay rates of the 
lightest Xi components (through higher dimensional operators) are typically much smaller than 
H{T) at freeze-out for both weakly and strongly interacting y. We thus have effectively a two 
step process and we can treat freeze-out and decay separately. 

In case y is a color singlet, the cosmological relic abundance will generically be determined 
by its (co)annihilations into EW gauge bosons resulting in a cosmological density of ~ 
0.01. On the other hand, the hnal relic abundance of a colored multiplet is determined in two 
stages. At temperatures T ~ m^/30 the relic abundance is determined by perturbative QCD 
annihilations resulting in ~ 10“^. Then, y undergoes a second stage of annihilation after 
the QCD phase transition, further reducing its relic abundance to a value roughly three orders 
of magnitude smaller [45]. 

The y lifetimes determine at which cosmological epoch they will decay. Such decays will 
involve the creation of energetic SM particles, which can produce a variety of observable effects. 
First, the decays of heavier multiplet components into the lightest y* state always happen well 
before nucleo-synthesis and give a negligible entropy release. On the other hand, if the lightest 
Xi states can decay through d = 5 operators, their lifetimes are at least of the order (0.1 — 10^) s, 
and may thus affect the primordial generation of light nuclear elements [46]. For longer lifetimes 
of the order (10^^ —10^^) s, y* decays would create distortions in the thermalization of the cosmic 
microwave background (CMB) before recombination. Such distortions of the spectrum by the 
injection of high-energy photons into the plasma lead to strong constraints [47]. Decays of 
Xi after recombination can give rise to photons that free-stream to us, and are visible in the 

the one-loop level there are no contributions from SM fermions and the charge conjugation C transfor¬ 
mation is a symmetry of the gauge and scalar sectors. The selection rules for these decay channels follow from 
the presence of this symmetry. 
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diffuse gamma ray background [48]. Observations by Fermi LAT [49] limit the flux of these 
gamma rays and thus constrain such scenarios. In general these observations of the diffuse 
gamma ray background rule out Xi with lifetimes between (10^^ — 10^®) s. If Xi only decay 
through d > 6 operators, they will survive to the present day. In case they are integer charged, 
they would act as heavy positively charged nucleons, producing anomalously heavy isotopes. 
A combination of measurements places severe limits on the abundance of terrestrial heavy 
elements today [50], effectively exclnding such scenarios. Alternatively, if their charge is a 
non-integer fraction of that of the electron, they are exclnded by the nnll results of searches for 
fractionally charged particles in bnlk matter on Earth or meteoritic material [52, 53]. Finally, 
snfhciently stable neutral Xi can form (a fraction of) dark matter, a possibility, which has been 
thoronghly covered in the literature [4-7]. In principle, one could think abont introdncing 
an inhnitesimal hypercharge (ey) which would make y absolutely stable but still pass all the 
cosmological bounds. This would open up additional DM candidate scenarios like the complex 
scalars (1,1, ey), (l,3,ey), (l,5,ey) or the Dirac fermions (1,1, ey), (l,3,ey). Representations 
having Y ^ 0 for ey —)■ 0 are excluded by direct DM searches [4]. On the other hand, higher¬ 
dimensional 811 ( 2 ) 2 , representations have a Landan pole below 10^® GeV (cf. Table 9). We will 
not entertain snch a possibility any fnrther since it is a rather simple distortion of the minimal 
DM setnp (see for instance [7]). We also refer the reader to existing literatnre for more details 
on experimental bonnds on ey (e.g. [54]). 

We close this section with a few general comments abont the possible interpretation of 
cosmological DM within our framework. First of all, we note that the microscopic natnre of the 
DM is still nncertain. For example Massive Astrophysical Compact Halo Object (MACHO) 
made of ordinary baryons (like black holes or neutron stars) could in principle be a viable 
option. It is known, however, that in snch cases departnres from the standard Big Bang theory 
are needed. The present cosmological data and various theoretical considerations favor the 
hypothesis of particle DM. Besides the minimal DM cases, repairing DM of this type in our 
framework means departing from minimality. The easiest possibility then is to assume, on top 
of the (non DM) accidental matter state, the presence of the fermionic minimal DM mnltiplet 
at ~ 10 TeV. However, this works only for some accidental matter states. In other cases extra 
d = 5 operators can trigger too fast decay of the minimal DM candidate. Additional possibilities 
inclnde the presence of extra gange interactions where the stability of DM is again guaranteed 
by an accidental symmetry of the new gange sector (for a recent work along these lines see [55]) 
or axion DM with PQ symmetry breaking above Aefr [56-58]. 

4.1 Relic abundance 

We hrst consider scenarios with uncolored y, where its lightest component is electrically charged. 
In these cases, direct searches already limit 3> mz (see Sect. 5.2) and we can compnte the 
relevant thermally-averaged cross-sections in the SU(2)i-symmetric limit. This approach is 
valid as long as all SU(2)i mnltiplet components are present in the thermal plasma. 

For inter-multiplet splittings of typical radiative size all heavier Xi components decay into 
the lightest one with lifetimes (cf. Fig. 5) which can be comparable or even shorter than the 
inverse Hubble rate at freeze-out (typically 0(10“^^ s)). Thus the abundance of the lightest Xi 
component (before itself starts decaying) is actually described by the sum of the densities of all 

note however that in principle these bounds can be evaded for »TeV [51]. 
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Xi states. And as long as Xi ^ Xj conversion rates are in eqnilibrium at freeze-ont (which is 
always the case for color singlet weakly interacting x), the actnal Xi Xj^ rates do not affect 
the total relic abnndance [59], and the SU(2)2, symmetric approximation can be justihed. 

Finally, we also ignore thermal corrections. They mainly indnce thermal mass splittings of 
the order Am^ ~ {g'^TY /which can be neglected at the level of precision we are considering 
here [5]. 

Due to the above approximations we can write a single Boltzmann equation that describes 
the evolution of the total abundance of all components Xi of fhe multiplet as a whole. In 
particular, it includes all co-annihilations in the form of 'Yhij ^AiXiXj SM particles). The 
hnal X abundance can be well approximated as [4] 


y ^ _ / fSO 1 ~ /oQl 

s(T) ^ Vvr^SMmpiT^(an)’ ^ ^ ^ 

where is the number of degrees of freedom of a whole x multiplet including anti-particles 
in case of complex representations, is the number of SM degrees-of-freedom in thermal 
equilibrium at the freeze-out temperature Tf (c.f. [60]), and s is their total entropy. The typical 
freeze-out temperature is Tf ~ m^/26 m^, such that we can keep only the dominant s-wave 

(co)annihilation processes. The relevant formulae for the corresponding thermally averaged an¬ 
nihilation cross-sections {av) into SU(2)^ 0 U(l)y vector bosons for both scalar and fermionic 
X with generic SU(2)^ 0 U(l)y quantum numbers can be found in [4]. The resulting es¬ 
timates are within 10% of the more complete treatment including p—wave annihilations and 
renormalization of the SM gauge couplings [5]. However, for rriy^ > 1 TeV, the relic abundance 
is expected to be further reduced by 0{1) non-perturbative (Sommerfeld) corrections due to 
the electrostatic Coulomb force effects [61]. In case of scalars, additional renormalizable Higgs 
portal interactions can also contribute to the annihilation cross section deferring freeze-out. In 
light of this our estimates of using Eq. (29) with dominant EW gauge boson contributions 
to {av) can be considered as upper bounds on the actual relic abundances of x- 

In the case of colored X; oue needs to consider two separate regimes of annihilation. The 
hrst era is before the QCD phase transition when all x components are freely propagating in 
the QCD plasma and the annihilation cross-section can be determined using perturbative QCD. 
The second era is after the QCD phase transition when the heavier multiplet components have 
decayed to the lightest X) which in turn have become conhned in color neutral bound states. 
The annihilation cross section in this second period turns out to be much higher than in the 
hrst, thus leading to a second period of annihilation which completely determines the hnal 
X relic abundance [45]. In particular, heavy colored particles are conhned within hadronic 
states of typical size i?had ~ GeV“^ which annihilate with a geometrical cross section yielding 
{av) ~ ^TR\^^^JTB/ra^, where Tb ~ 180 MeV is the temperature at which QCD conhnes and 
hadronic bound states form. The hnal x relic abundance can thus be approximated as 
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GeV-^ 
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(30) 


where we have used Eq. (29) with Tf = Tb and gsM ~ 15 just below the QGD phase tran¬ 
sition. The annihilation proceeds through intermediate excited bound states which decay by 
radiating away photons before annihilating into quarks and gluons [45]. These processes need 
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Figure 10: Comparison of NDA estimates for the Xo lifetimes (r^) in case of colored multiplets 
(8,1, 0)5 (drawn in thick dashed purple), (8, 3, 0)5 (drawn in thin brown) and (8, 3,1)5 (drawn 
in thick green) with the annihilation lifetimes of the corresponding xo hadronic bound states 
('Thado, drawn in thin dashed gray). In the shaded region, the xo decay before their hadronic 
bound states fully annihilate. See text for details. 


to be considered carefully, since such late decays to photons and hadronic jets could affect 
nucleosynthesis [62]. In case of electrically charged Xi, this process is fast with a lifetime of 


'hiad+ ~ 


ttEMAhad 


~ 3 X 10 “^^ s 


a^(mx) / Ahad A ^ 

0.1 J VGeVy VTeV/ 


(31) 


where Ahad ~ 1 GeV is related to the QCD string tension a via a ~ A^^^. On the other hand, 
for electrically neutral y*, radiation of photons is loop suppressed, leading to a much longer 
annihilation process 


47rm^ / Tb\ x 6 \ 

“InALd Uhad; ^ n2.7TeVy' \GeV) ^80 MeVy* 


(32) 


where the inequality is due to neglected non-local contributions to the decay rate. This scenario 
however, only applies to our cases x ~ (8,1, 0)5, (8, 3, 0)5 and (8,3, 1)5 when the lightest 
component is neutral. In Fig. 10 we plot the relevant lifetimes and Thado as a function 
of X mass. We observe that in the low mass regime, Thado < t'x and we can use the non- 
perturbative result in Eq. (30) to estimate the hnal x abundance. However in the heavy x 
limit, ThadO > W and so x decays before the second stage annihilation process is completed. 
In this case the relevant abundancies are those after the hrst stage of annihilation above the 
QCD phase transition. We can estimate them using the general SU(iV) annihilation cross- 
sections computed in [63] after exchanging the relevant SU(iV) group invariants and correcting 
for the different number of degrees of freedom. In particular for the adjoint representation of 
QCD Y.i,j,a,b = 216 and I[D“, D^]jiP = 72. Velocity expanding the resulting 
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XX 99 cross section we obtain (c.f. [64]) , 


{av) 


277ra‘l 

9^m? ’ 


(33) 


in the conventions of [4] . In estimating the resulting relic abundance we can safely assume that 
weak interactions keep Xi ^ Xj processes in equilibrium until decoupling [59]. We have also 
checked that employing Eqs. (29) and (33) and using the running as{2m^) to estimate the relic 
abundance reproduces the results of the full leading order perturbative QCD calculation and 
integration of the Boltzmann equation to within 20% (in agreement with similar results for 
the fundamental QCD representation in [63]). Non-perturbative Sommerfeld corrections are 
expected to lead to an 0{1) reduction in the hnal result and so our estimates can again be 
taken as upper bounds on the actual relic abundance of color octet scalar y above the QCD 
phase transition. In the intermediate regime when Thado W the actual evolution of the 
X number density in the primordial plasma depends on the detailed dynamics of the neutral y 
hadron annihilation and decays, the evaluation of which is beyond the scope of our study. As 
we show in the next section however, these details are never relevant, since they do not lead to 
observational constraints. 


4.2 Implications for Big Bang Nucleosynthesis 

In general, nucleosynthesis of primordial elements in the early Universe represents a sensitive 
probe of any metastable relic with lifetime of about 1 s or longer [46, 65]. The constraints come 
from two classes of processes: injection of very energetic photons or hadrons from decays during 
or after BBN adds an additional non-thermal component to the plasma and can modify the 
abundances of the light elements [66-70]; in addition, if the relic particle is electromagnetically 
charged, bound states with nuclei may arise that strongly enhance some of the nuclear rates and 
allow for catalysed production of e.g. ®Li, ^Li [71-73]. The Standard BBN prediction for the 
®Li abundance is actually signihcantly smaller than the observed one, so that the presence of a 
charged relic with appropriate lifetime can help reconciling BBN with the measured abundances 
of ®Li and "^Li [74-79]. 

In general, the decay can produce very energetic SM particles that can initiate either 
hadronic or electromagnetic showers in the plasma. The most stringent bounds are obtained for 
a relic that produces mostly hadronic showers, since electromagnetic particles like photons or 
electrons can thermalize very quickly by interacting with the tail of the CMB distribution until 
times of about 10® s. In the following we will consider the constraints for relics producing a 

small number of energetic hadronic jets with a branching ratio Shad = 1 and Shad ~ where 

Shad is the decay energy released in the form of hadronic showers to obtain conservative upper 
bounds on y number densities. This assumption is mostly valid if x is colored (in particular 
in this case always Shad = 1 and Shad ^ 'nr^/2), while Shad < 1 is expected for non-colored x- 
Then the hadronic BBN bounds are relaxed accordingly by a factor 1/Shad- Finally, for life¬ 
times > 10^ s, electromagnetic interactions start having a signihcant effect and the bounds 

above > 10^ s become effectively independent of Shad- 

In practice there are three regions of lifetimes as discussed in [62]: for lifetimes 0.1 s < < 

100 s the dominant effect is the interconversion between protons and neutrons, that changes the 
^He abundance by overproducing it; at longer lifetimes 100 s < < 10^ s hadrodissociation 
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is the most efficient process and the bounds come from the non-thermal production of lithium 
and deuterium; hnally at late times 10^ s < < 10^^ s photodissociation caused both by direct 

electromagnetic showers and by those generated by the daughter hadrons starts to dominate 
and result mainly in the overproduction of ^He. 

In the following we use the results from the general analysis of [62] for the second and the 
third lifetime regions. In particular the bound coming from the abundance of ^He scales as 
1/m^ and does not depend on ^Bhad (we neglect the decay energy released into neutrinos, which 
is always expected to be a small fraction of m^). On the other hand due to the Li anomaly, 
in the second region we only consider bounds on coming from the deuterium to hydrogen 
abundance ratio (D/H), which scale roughly as Finally, we note that a charged thermal 

relic with ~ 10^ — 10^ s and abundance just below the D/H bound may (partly) ameliorate 
the standard BBN Lithium problems (c.f. [80]). 

The value of the observed ^He abundance Yp = 4(nHe/’^H)/(l + duHe/'^.n) ~ 2(?7,„/np)/(l + 
nn/^p) which dominates the constraints in the hrst lifetime region has been updated since the 
analysis of [62] and currently reads Yp = 0.250(3) [65] to be compared with the prediction of 
standard BBN of I^sbbn _ 0 2483(5) [81]. The bounds from [62] which assumed signihcantly 
smaller Yp thus need to be re-evaluated. For this purpose we numerically solve the relevant 
Boltzmann equations (c.f. [60]) 
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^pnVp ^npVn , 
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where i/i = rii/nb and Ub is the baryon number density (we use rj = Ub/n^ = 6.1 x 10 ^°), 
H{T) = /m■p\)^/gsuJ^ is the Hubble rate and Tn = 880(1) s is the neutron lifetime. 

We have furthermore dehned = Kn^p — Kp^n and Xnp — Apriexp(l/?/) ~ (1443/r„)?/^(?/ -|- 
0.25)^ [82, 83], where y = T/Q and Q = 1.293 MeV is the neutron-proton mass difference. 
We have checked that this approximate form of Xnp and Xpn reproduces the hnal results using 
exact numerical integration of the weak nucleon conversion rates (c.f. [82, 84]) to better than 
0.5%. Finally, for the catalyzed nuclear conversion rates we employ the formulae for Km^n' 
including all the numerical inputs as dehned in [62]. At temperatures much bigger than Q we 
expect yn = {1 — yp) = 1/(1 -|- exp(Q/T)) and y^ given by its thermal relic abundance y^. 
Finally, the resulting ^He abundance is well determined by yn^p at (T ~ 8.5 x 10® K) when 
BBN begins [81, 85]. As a cross-check of our approach, we have determined Yp in absence 
of y^ and obtained Yp{y^ = 0) = 0.243, which is consistent with the expected precision in 
light of our approximations. In particular, neglected higher order effects would increase Yp by 
2% [82], reproducing the standard BBN result. In setting constraints we use our estimates only 
to compute the deviations of Yp from the standard BBN value AYp = Yp{y^) — Yp{y^ = 0) and 
compare Ah/, to the 2a region of the observed Yp value. As discussed above, the effects 

of X decays on Yp scale as 1/Bhad, and for Bhad ^ 100 GeV (in the form of a hxed number of 
hadronic jets) also approximately as 

We hnally determine the upper bound on the possible contributions of x decays to Yp 
by hxing Bhad = 1 and assuming x decay to two hadronic jets. The resulting ehects then 






Figure 11: Comparison of BBN Yp (left hand panel) and D/H (right hand panel) constraints 
on the abundances and lifetimes of metastable hadronically decaying particles with the cor¬ 
responding estimates for the cases of viable x multiplets. Each line, going from bottom to 
top, corresponds to the mass range 0.5 TeV < < 5 TeV. The only explicitly labeled 

examples are (1, 5, 2 )s,Qlp= 3, (1, 5, 2)s,Qpp=4, (1, 5, l)s,Qpp=3 and (1, 2, 5/2)5,Qpp=3, which are 
potentially constrained. The case (1, 7, 0) is not shown as its lifetime is longer than 10® s for 
0.5 TeV < < 5 TeV and Aeg = 10^® GeV. 

scale as ' and we can use a single reference value to constrain the abundances of y at 
different m^. As discussed above, this leads to a conservative 0{1) overestimate of the actual 
Yp constraints for decays of x involving also uncolored hnal states. The comparison of the 
lifetimes and relic abundances of all x candidates from Tables 1 and 2 in the mass range of 
0.5 TeV < < 5 TeV with the Yp bound estimated in this way is shown in Fig. 11 (left hand 

side). Notice the almost discontinuous drop of the abundance for some representations. This is 
due to the fact that for colored multiplets featuring a neutral component there is a qualitative 
change of behaviour when the lifetime of the particle becames smaller than the annihilation 
lifetime of the associated hadronic bound state (cf. Fig. 10). We observe that all the x are 
consistent with the Yp constraint. Also, most of the candidates have lifetimes shorter than 
~ 10 s (for Agff ~ 10^® GeV), so that no further bounds from BBN processes at later times can 
be derived. The only exceptions are the cases in Table 6 where the lightest x component can 
only decay through long cascades involving off-shell heavier components and W bosons as well as 
(1, 7, 0)5 decaying exclusively through loop-induced processes. For these cases the D/H bound 
applies as shown in Fig. 11 (right hand side). In particular, while all the colored multiplets 
(including those decaying with cascades) are consistent with this constraint due to their low 
relic abundance after the second stage of strong annihilations, all the long-lived uncolored cases 
are in general constrained. In the relevant region of relic abundances and lifetimes, the bound 
turns out to be insensitive to the exact x relic abundance or decay mode and so even our crude 
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Figure 12: Dependence of the D/H bound on the x mass as a function of the EFT cut-off 
scale Aeff suppressing the relevant d = 5 decay mediating operators. See text for details. 


estimates suffice to extract fairly robust lower bounds on y masses. They are shown in Table 7 
for a hxed value of Aeg = 10^^ GeV, while the Aeg dependence is shown explicitly in Fig. 12. 

We hnally note that it is close to these D/H exclusion bounds where the primordial Lithium 
problem might be addressed by the presence of y listed in the hrst three rows of Table 7. 
A detailed exploration of this possibility goes however well beyond the scope of the present 
analysis and we leave it for future study. 


Spin 

A 

Qlp 

Mass bound [GeV] 

0 

(1,2, 5/2) 

3 

790 

0 

(1,5,1) 

3 

920 

0 

(1,5,2) 

3,4 

530, 1900 

0 

(1,7,0) 

0,1, 2, 3 

> 5000 


Table 7: BBN bounds on the masses of long lived y multiplet components assuming hxed 
Aeff = 10^5 GeV. 


5 Collider phenomenology 

In this section, we discuss the collider phenomenology of the new states and give bounds on 
the masses of the new particles. The states of Table 1 can be ordered into two classes: the ones 
which can decay by renormalizable interactions and the ones which decay via effective operators. 
In case the states decay via renormalizable interactions, they can be detected via their decay 
products. We will shortly comment on the renormalizable cases in Sect. 5.1, mostly referring 
to the existing literature. If they only decay via effective operators they are rather long-lived 
and can eventually leave the detector before decaying. If the new particles are uncolored, the 
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signature depends on whether the particle is charged or not. A summary of the different mass 
bounds for the uncolored cases can be found in Table 8. 

All the extra colored states, given in Table 2, can only decay via d = 5 operators and are 
hence long-lived. They hadronize and build exotic new mesons or baryons. We will discuss 
them in Sect. 5.4. 

The production of the new exotic fermions and scalars proceeds via Drell-Yan processes. 
Throughout this section, we use the LO Drell-Yan production cross sections. Formulae are given 
e.g. in Ref. [37]. The cross section for scalars is in general more than one order of magnitude 
smaller than for fermions, which explains the lower exclusion bounds on the scalars. 


Spin 

X 

Qlp 

Mass bound [GeV] 

0 

(1,2,3/2) 

1,2 

430, 420 

0 

(1,2,5/2) 

2,3 

460, 460 

0 

(1, 5, 0) 

0, 1, 2 

75, 500, 600 

0 

(1, 5, 1) 

-1, 0, 1, 2, 3 

640, 50* (85), 320, 490, 600 

0 

(1, 5, 2) 

CO 

r—1 

o' 

85, 530, 410, 500, 570 

0 

(1, 7, 0) 

0, 1, 2, 3 

75, 500, 600, 670 

1/2 

(1.4,1/2) 

-1 

860 

1/2 

(1,4,3/2) 

0 

90 

1/2 

(1, 5, 0) 

0 

95 


Table 8: LHC-I/LEP summary bounds for uncolored accidental matter multiplets decaying 
via d > 5 operators. Bounds on the neutral particles are given under the assumption of very 
small mass splitting to the |Q| = 1 component. The exclusion bound in braces corresponds to 
the case where the next-to-LP has Q = —1 instead of Q = 1. *A stronger exclusion bound, 
depending on the size of the portal coupling a (see Eq. (6)), can be obtained from the Higgs 
data. 


5.1 Renormalizable cases 

Among the extra multiplets which preserve the flavor group of the SM, compatibly with cos¬ 
mology and a cut-off scale of Aeg ~ 10^® GeV, we identihed four states which decay via renor- 
malizable interactions, namely (1,1, 0)5, (1,3,0)5, (1,4,1/2)5 and (1,4, 3/2)5. In all of these 
cases, the new state acquires a VEV. With the only exception of the SM singlet, these VEVs 
must be small in order to comply with EW precision measurements. For 0{1) couplings in the 
scalar potential and barring hne-tunings this implies >2 — 20 TeV, cf. Sect. 2.2.1. 

The gauge singlet can sizeably mix with the Higgs boson. Such mixing is constrained by 
the current Higgs data, see e.g. Ref. [86]. The triplet and quadruplet scalar multiplets can only 
have a very small mixing with the Higgs boson due to its effects on EW precision observables. 
Nevertheless, their charged components can modify the Higgs to 77 and to Zj rates by their 
loop contributions. Whether these loop contributions suppress or enhance the diphoton rates 
depends on the sign of the couplings in the scalar potential [87, 88]. Finally, masses of the new 
neutral scalars below 62.5 GeV can be probed by the invisible Higgs boson width. For more 
details see Sect. 2.2.2. 

Low masses of the triplets and quadruplets can be constrained by the Z width (see Sect. 5.3 
for more details). Their charged components can also be directly detected. They decay into 
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vector bosons or via cascades into vector bosons and the (off-shell) neutral components of the 
multiplet. The coupling to two vector bosons is proportional to the VEV of the multiplet. Apart 
from searches for singly charged Higgs bosons, searches for multiple charged Higgs bosons can 
provide a distinctive probe for large scalar multiplets. By now, searches for doubly charged 
Higgs bosons [89, 90] have only been performed for decays of the charged scalars into fermions, 
as in e.g. the case for models with Y = 1 triplets [91] and quadruplets with additional vector- 
like matter for seesaw mass generation of neutrinos [92, 93]. Bounds on doubly charged Higgs 
bosons decaying to can be obtained by reinterpreting SUSY searches for dileptons, 

missing energy and jets, and can exclude masses of the doubly charged scalars up to roughly 
200 GeV at = 7 TeV for SU(2)i triplets [94], 

5.2 Colorless and charged LP 

Charged stable particles will undergo charge exchange with the detector material. For masses 
larger than 100 GeV the time of flight till to the outer detector is signihcantly larger than 
for lighter objects such as muons. They can hence be distinguished by their longer time of 
flight and by their anomalous energy loss in the detector. The energy loss in the detector is 
described by the Bethe-Bloch formula and depends on the speed and the charge of the particle. 
Searches for such ionizing tracks have been performed in Refs. [95-108]. The strongest bounds 
come from the CMS search of Ref. [108], where the exclusion limits on the production cross 
sections of fractionally, singly and multiply charged particles are presented assuming vanishing 
quantum numbers under SU(2)Hence, in order to use the results of Ref. [108] they need to be 
recast. Reference [109] gives tabulated efficiency values in terms of the transverse momentum, 
the pseudo rapidity and the velocity (3 of the heavy charged particle. These can be used to 
reinterpret the results of Ref. [108] without running a full detector simulation. 

In order to compute the efficiencies, the models were implemented into Madgraph 5 [41] with 
the help of FeynRules [42]. The cross sections, computed at LO at the scale Q = \/§ using the 
MSTW2008 [110] parton distribution functions, were than rescaled by a factor accounting for 
the change in the efficiencies with respect to the cases considered in Ref. [108]. 

We hud that the efficiencies for fermions with non-vanishing SU(2)£ quantum numbers 
barely change compared to the case with T3 = 0. The exclusion limits from Ref. [108] on 
the cross section can hence be applied naively. For scalars, the efficiencies change slightly 
compared to the fermions. For masses of 300 GeV the efficiency is slightly smaller than the one 
for fermions, for 800 GeV it is roughly 15% larger. For the cases we considered, the efficiency 
for scalars is always within 3% of the case with vanishing quantum SU(2)j;, quantum numbers. 

In order to derive exclusion limits for the scalars, we adopt the following procedure. We 
compute the cross section and derive an approximate bound using the 95% G.L. upper limits 
given in Ref. [108]. With this approximate bound at hand, we compute at the naive bound the 
efficiencies for the scalar and compare it to the efficiencies of a fermion with T3 = 0. The results 
can then be recast by the appropriate factor. We note, however, that such rehned bounds are 
in good agreement with the results obtained naively. 

Fven though the efficiency values in Ref. [109] are given for the singly charged analysis only, 
we apply the same procedure to the |Q| > 1 case to check whether also here, the naive method 
gives sensible results, as the basic cuts in both analysis are the same. Indeed we hnd also here 
that within the precision of our results, the naive estimate is very good. The results can be 
found in Table 8. The limits for charged fermions are stronger than for scalars due to the larger 
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production cross section. The weakest exclusion limits are obtained for T3 = 0. 

Before concluding this subsection let us mention the recently approved LHC experiment 
MoEDAL [111, 112], whose target is the study of new physics phenomena (e.g. magnetic 
monopoles) which manifest themselves through the presence of highly-ionizing particles. In 
particular, the nuclear track detectors of MoEDAL are sensitive to particles with \Q\/(3 > 5, 
where Q is the charge and (3 is the velocity of the particle in units of the speed of light. For our 
framework, with Q ranging from 1 to 4 (cf. Table 8), the discovery potential of the MoEDAL 
experiment will be relevant at low values of the f3 distribution. 

5.3 Colorless and neutral LP 

The search for stable (on detector scale) neutral and colorless particles is very challenging at 
the LHC. Limits can either be set directly on the mass of the neutral particle as e.g. by mono-x 
searches or from constraints on the invisible Z width or, indirectly, by giving bounds on the 
mass of the second lightest particle of the multiplet, such as in disappearing track signatures. 
Let us discuss in turn all these possibilities: 

• Mono-x searches 

Neutral stable particles are searched for at the LHC in mono-x searches, in which large 
missing energy is accompanied by a radiation of an additional high-energetic particle x, 
where “x” can stand for a jet, a photon, aW 01 Z boson, a top quark or a Higgs bosons. 
Nevertheless, we hnd that the monojet searches of Ref. [113], which potentially have the 
strongest reach [114], are not sensitive to our states yet. Similar results were found for 
instance in Ref. [7], in the case of a fermionic (1,3,0) multiplet. Monojet searches can, 
however, become sensitive at 14 TeV [7]. 

• Invisible Z width 

At LEP, the Z boson width was determined with high accuracy [115, 116]. This measure¬ 
ments set a tight bound on new physics contributions to the invisible Z width at the level 
of T?]))^ < 2 MeV. This hence excludes charged particles, or particles with non-trivial 
SU(2)i quantum numbers, up to the kinematic bound for the Z —>■ XX^ decay, meaning 
that masses m-^ < 45 GeV are excluded. 

• Disappearing tracks 

Disappearing tracks can be observed at the LHC if a rather long-lived charged particle 
decays within the sensitive volume into a neutral particle and a soft pion, which is not 
detected. The strongest limits on these searches [117], are sensitive to lifetimes of the 
charged particle between 0.1 ns and 10 ns. We checked whether the typical lifetimes for 
our particles lie within this range. It turns out, however, that for all fermionic states of 
Table 1 with a lightest neutral state, the mass splitting between the neutral component 
and the charged component is always so large, that the lifetime is smaller than 0.1 ns 
(cf. Fig. 5). This is due to the fact that the radiative mass splitting increases with the 
hypercharge and the SU(2)i quantum number. For the scalar states featuring a lightest 
neutral component and which do not decay through renormalizable interactions, the same 
argument holds. In addition, the mass splitting does not need to be purely radiative but 
a larger mass splitting can also stem from the potential term. 
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• LEP bounds on charginos 

The LEP experiments set bounds on charginos that are nearly mass degenerate with the 
lightest neutralino. These bounds can be reinterpreted for our purposes in order to derive 
limits on the mass of the lightest neutral particles, since we showed in Sect. 2.3 that 
the next-to-LP has always charge Qlp i 1- References [118-121] cover a mass splitting 
Am between 200 MeV < Am < 5 GeV and are based on soft events with an initial 
state radiated photon. In order to estimate the limits for the case where the LP of the 
multiplet is neutral, we took the OPAL results of Ref. [121]. There, the results were given 
in terms of a 95% C.L. upper limit on the cross section. We implemented the models 
into MadGraph 5 with the help of FeynRules [42], computed the cross section values for 
the IQI = 1 charged component of the multiplet, and compared them to the given limits 
in Ref. [121]. In order to verify this procedure, we computed the efficiencies for example 
points on parton level, using MadAnalysis [122]. For the fermionic states we found that 
the efficiencies are basically unchanged compared to the chargino case. For the scalars 
they turned out to be a bit reduced, which is however not relevant given the precision 
to which we estimate the bounds. The hence obtained limits on the charged components 
are given in Table 8 . At the accuracy we are working this essentially corresponds to the 
bounds on the neutral components, which are obtained after subtracting the small mass 
splitting. 

For the case of the scalar (1,5,1) multiplet, either the +1 charged or the —1 charged 
component can be the next-to-LP. These two cases lead to different exclusion bounds. 
The exclusion bound from the Q = 1 state being the second-lightest component is much 
smaller due to the smaller production cross section, as this state corresponds to = 0 . 
In such a case a stronger bound can come from the Higgs invisible width. 

What about mass splittings larger than 5 GeV not covered by the chargino search of 
Ref. [121]? If the neutral LP of a scalar multiplet is the component with the small¬ 
est/largest isospin, then the mass splitting between the IQI = 1 next-to-LP and the 
neutral LP can also be larger than 5 GeV. This is not true for a generic value of the 
isospin —j<I<j, as for the neutral state to be the lightest a cancellation between the 
tree-level and radiative mass splitting is required. In particular, the only case in which 
we have to consider a mass splitting larger than 5 GeV is for (1,5, 2 ) 5 . Even in such a 
case, however, mass splittings larger than about 20 GeV are excluded by EW precision 
observables (cf. Sect. 2 . 2 . 1 ). 

Searches for charginos decaying into neutralinos and W bosons (with the Ws decaying 
hadronically, semileptonically or leptonically) and for Am > 5 GeV were performed also 
at LEP [123].^^ In the (1,5, 2)5 case, for mass splittings around 5 GeV the estimated 
bounds on the |Q| = 1 particle decaying to the neutral state turn out to be weaker (by 
roughly 10 GeV) than the ones given in Table 8 . For mass splittings larger than 15 GeV 
the limits become tighter (by around 5 GeV). A more detailed analysis is however beyond 
the scope of this paper. 

^^LHC searches for pair production of charginos, with the charginos decaying to IPs and neutralinos, are not 
yet sensitive to such low mass splittings. 
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5.4 Colored LP 


The description of long-lived colored particles is complicated by the effect of non-perturbative 
QCD interactions. In fact, in all the cases of Table 2 the decay of the new states is induced by a 
d = 5 operator and the lifetimes are long enough that heavy colored particles hadronize before 
decaying. The theoretical description of the hadron formation and of the nuclear interactions of 
such states with matter represents the main source of uncertainty. In this section, we will briefly 
describe the various steps for the path of our new states, from their production to their escape 
from the detector or, in case the initial velocity jd is small enough, to their eventual stopping 
and decay inside the detector. We finally conclude by looking at the recent LHC results that 
are tuned to the case of QCD bound states of SUSY particles with quarks and gluons. We refer 
to [124] for a review on various phenomenological aspects of stable massive particle at colliders. 

Production. The production mechanism of the new particles in our framework is determined 
by the color quantum number and by the value of the mass. We are interested in the 
fundamental, the two indices symmetric and the adjoint representations of SU(3)c and we 
denote such cases respectively as Cs, Cq and Cg. At the renormalizable level the presence 
of a U(l) or a Z 2 accidental symmetry guarantees these states to be pair produced. 
The fate of the produced state crucially depends on the velocity at the production time. 
Relativistic particles will lose energy throughout the detector but eventually escape it, 
while slow particles will be stopped in the detector and decay at a later time. Typical 
velocity distributions are displayed in Fig. 13. There is no significant difference between 
scalars and fermions. For higher center-of-mass energies and lower masses, m^, higher 
velocities are more probable. Notice also that due to the normalization with respect to the 
total cross section there is no difference between the different color quantum numbers. 
Another issue that has to be considered for a complete description of the production 
mechanism is the Sommerfeld enhancement due to ladder exchange of gluons. This effect 
is relevant only for slowly produced states. 

Hadronization. Once produced, a colored particle combines with quarks and gluons to form 
a colorless hadron state. For example, a color-triplet can form bound states such as C^q 
or C 3 qiq 2 , an octet state can form invariants such as Csqq, Csqiq 2 q 2 or Cgg, while the 
sextuplet can hadronize in states of the form C^qg, Ceqqq and C^qq. The probability 
for Ci of hadronizing in a given bound state are rather uncertain, different models give 
quite different values. Bound states made of larger numbers of quarks and gluons are 
expected to be heavier [128] and, even if the hadronization in this channels could be non- 
negligible, the newly produced state could decay into a lighter one plus ordinary baryons 
and mesons, through QCD gauge interactions. Despite the fact that the hadronization 
processes are very uncertain, in some cases a detailed knowledge is not very important. 
As we are going to comment soon, nuclear conversions can wash out the information on 
the original hadron state at production. 

Propagation through matter. As soon as produced and during their propagation, the long 
lived colored particles interact with the electrons, protons and neutrons present in the 

^^See for example the comparison between the string model adopted by Pythia [125] and the cluster model 
used in HERWIG [126, 127] in Table 3 of [124]. 
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Velocity distribution colored fermion at Vi=8 TeV 


Velocity distribution colored fermion at Vi=i4 TeV 




Velocity distribution colored scalar at Vr=8 TeV 


Velocity distribution colored scalar at Vi=14 TeV 
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Figure 13: Velocity distribution normalized to the total cross section for colored fermonic 
states (top) and colored scalars (bottom). Left (right) panels correspond to a/s = 8 TeV 
(a/s = 14 TeV) for = 500 GeV (blue), = 800 GeV (pink) and = 1.1 TeV (light 
blue). 


detector. The dominant interactions we consider here are the electromagnetic and the 
strong ones. 

• Electromagnetic interactions. 

A particle with electric charge can interact with atomic electrons as well as with 
protons and neutrons in the nuclei. In the hrst case the net effect is the ionization 
of the atoms while interactions with atomic nuclei generate displacements of atoms 
from the lattice. In both cases the heavy long lived particle loses energy, however 
the energy loss dE/dx from ionization is much larger than the one generated in the 
non-ionizing way. 

The main role of the heavy elementary parton is to contribute to the electric charge 
of the hadron. Indeed the SM gauge quantum number and the requirement to obtain 
a color-singlet hadron state has an important influence on the total charge of the re- 


36 





















suiting bound state. We notice that considering all the cases in Table 2, the resulting 
bound state has always integer charge. The most unfavourable situation, from the 
point of view of detection, happens when the resulting hadron is electromagnetically 
neutral. 

• Strong interactions. 

Before discussing the various kind of interactions between the heavy hadrons and 
the matter in the detector, let us clarify the role of the parton Ci in the nuclear 
reactions. Due to their large mass, the wave-functions of the CiS are expected to be 
highly localized as compared to those of the light constituents (quarks and gluons) 
that are spread in space as in ordinary QCD. From this observation we can draw 
the conclusions that the probability for the heavy parton Ci to interact with matter 
is very low, while the typical cross section of the hadron with matter, being due to 
the effect of the light partons, is expected to be of the same order of those for pion 
scatterings. 

Our heavy long-lived hadron can have elastic as well as inelastic reactions with 
nucleons. Elastic scatterings are not particularly relevant, indeed the energy loss 
is small because the long-lived hadron scatters on a much lighter target nucleus. 
Inelastic processes are instead those responsible for the slowing-down of the hadron. 
In an inelastic reaction it is also possible to exchange baryon and electric charge. 
The importance of baryon exchange has been emphasized in [129], in these reactions a 
heavy hadron is transformed into another one with different baryon number. For the 
case of the gluino it has been argued that, bound states with null baryon charge {R- 
mesons) are very efficiently converted into baryonic states in reactions like {Csdd) -|- 
p ^ 71 + (Csudd). The reverse reaction is suppressed mainly because of the mass 
split ordering of the various hadrons and by the low presence of pions as targetd in 
the detector material. As a consequence, early in the detector, mesons are converted 
into baryons. 

Processes with charge exchange are particularly relevant for detecting the presence 
of the heavy parton. Indeed, tracks generated by the passage of electric particles 
can be easily detected, so it is important to understand the value of the electric 
charge of the hadron through all its travel in the calorimeter material. Even in the 
most pessimistic case of a hadron generated as a neutral bound state, reactions with 
charge exchange can covert it in a charged state that could be detected. 

We hnally mention that there are a series of phenomenological approaches to describe 
the strong interactions between heavy hadrons and matter [130-132]. Despite the 
fact they agree on several general qualitative aspects, they give rise to rather different 
quantitative results. 

Stopping and decay. Depending on the Ci mass, a non-negligible number of particles could 
stop in the detector. In our framework most of the states are supposed to decay within 
few seconds. This constitutes a really interesting possibility to understand the structure 
of the d = 5 effective operators responsible for the decay. However, the detection of these 
processes represents a severe experimental challenge. 

Having described the most important aspects of the phenomenology of our long lived par¬ 
ticles, we now move to comment about direct searches. At the LHC, searches for long-lived 
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colored particles are performed in the context of -R-hadrons, which are bound states of gluinos/ 
squarks and quarks/gluons. The i?-hadrons can be detected by the longer time-of-flight to the 
outer detectors and their anomalous energy loss. Complementary to the searches relying simply 
on the longer time-of-flight and the anomalous energy loss, are searches for stopped i?-hadrons. 
They are particularly suited for velocities (3 <^1. 

In the former case, the strongest limits come from the CMS search of Ref. [108], which 
excludes gluino masses up to 1276 GeV, if the fraction of gluinos hadronizing into g — g bound 
states is 0.5. If the fraction is equal to one, gluino masses are excluded up to 1250 GeV. In 
such a case the R-hadron is neutral in the inner tracker which leads to a smaller energy loss 
and hence a lower exclusion bound. In Ref. [108] stop masses were excluded up to 935 GeV 
(818 GeV). The exact exclusion bound depends on the modeling of the interactions of the stop 
with the detector material. The exclusion bound is here given for the so-called cloud model 
of Ref. [129, 133] (the charged-suppressed model of Ref. [134]). Very similar exclusion bounds 
were obtained in a recent ATLAS study [135] both for long-lived stops and gluinos. 

In the search for out-of-time decays of stopped gluinos or squarks of Ref. [136], gluinos 
are excluded up to masses of 880 GeV, assuming BR(^ —)■ x^g)=100% and a gluino lifetime 
between 1 /rs and 1000 s. Stop masses are excluded up to 470 GeV for BR(f —)■ x^t)=100% 
and stop lifetimes between 1 /rs and 1000 s. The exclusion bounds require that the neutralino 
mass is kinematically consistent with the used requirements on the energies of the gluon or 
respectively top decay products. The search furthermore assumes a cloud model for the R- 
hadron interactions. Reference [137] excludes sbottom masses up to 344 GeV for BR(6 —)■ 
x)’&)=100% for lifetimes between 1 /is and 1000 s. 

The results on these searches of R-hadrons cannot straightforwardly be applied for our cases. 
In order to give exclusion bounds on our states a detailed study of the hadronization of the 
different states and a full detector simulation would be necessary. We expect that the exclusion 
limits not only depend strongly on the SU(3)c quantum number, but they depend also on the 
charge and 811(2)^ quantum number. The charge influences the energy loss directly, whereas 
different SU( 2 ) 2 , quantum numbers lead to different time intervals, in which the R-hadron 
propagates as a neutral particle in the detector due to pion exchange between the different 
members of an isomultiplet [138]. For stopped R-hadrons, in addition, the searches depend on 
the BRs. For our cases if the new exotic particle decays into a missing energy and jet signature, 
there are always several other operators allowing the particle to decay (cf. Table 2), such that 
the corresponding BR likely deviates from 1. 

A further discussion on the bounds on long-lived colored states is beyond the scope of this 
paper, as, even in the well-studied case of R-hadrons the exclusion bounds depend signihcantly 
on the modelling of the hadronization and the nuclear scattering model. 

6 Conclusions 

Low-energy tests of fundamental symmetries provide a powerful probe of new physics scales up 
to energies of about 10^® GeV. Given the accidental (B and L) and approximate (GP, flavor 
and custodial) symmetry structure of the SM, it is somewhat surprising that signals of physics 
beyond the SM (if it exists) have not been observed so far. 

This last statement hinges on the theoretical prejudice that new physics effects at low 
energies can be described by a generic EFT, where the Wilson coefficients of the effective 
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operators are 0{1). There are, of course, exceptions to this point of view. The simplest one is 
maybe to allow for ultraweak couplings in the theory, so that the generic EFT power counting 
fails - an extreme example being the dissolution of the Weinberg operator when RH neutrinos 
are below the EW scale. On the other hand, the effective operators might not be there due to 
an exact symmetry of the Lagrangian in the full theory (as e.g. B — L in left-right symmetric 
models [139]) or they might be suppressed due to an approximate symmetry (as e.g. in minimal 
flavor violation [140]). In this paper, we explored yet another possibility: the quantum numbers 
of the new physics states below the EFT cut-off are such that by only requiring Lorentz and 
SM gauge invariance, the accidental and approximate symmetries of the SM are automatically 
preserved at the renormalizable level. The resulting new physics dynamics is practically invisible 
to low-energy indirect searches, and the only way of experimentally probing these scenarios is 
by direct production and detection of new particles at colliders. We hence focused on the 
phenomenological possibility that the new states lie within the kinematical reach of the LHC. 

Barring few exceptions, the new matter multiplets are subject to extra accidental Z 2 or 
U(l) symmetries which forbid their decays at the renormalizable level. Whenever the LP in 
the multiplet is color- and charge-neutral, it forms a DM candidate [4-7]. Generally however, 
the extra multiplets will decay due to the presence of higher dimensional operators in the EFT. 
In the spirit of generic EFT we choose to work with a cut-off scale of Aeg ~ 10^® GeV which 
is large enough not to require any further protection mechanisms in the full theory and is 
moreover suggested by the observations of neutrino masses. The inhnite set of possible states 
which satisfy the accidental symmetry conditions can then be reduced thanks to cosmological 
considerations. In particular, since the new states are long-lived, scenarios where the lightest 
component of the new multiplet is charged and/or colored are constrained by cosmological 
observations as well as by searches for exotic forms of matter on the Earth and in the Universe. 
The latter practically exclude all the cases where the charged and/or colored LP decays via 
d > 5 operators. 

Another handle in order to further reduce the list of possible states is the requirement that 
the theory remains perturbative up to the cut-off scale of the EFT. In particular, we required 
that no Landau poles are generated below Ae® ~ 10^® GeV. As a byproduct of the perturbativity 
analysis we noticed that, due to accidental cancellations in the coefficients of the one-loop beta 
functions for the non-abelian gauge factors, two-loop corrections can become important and 
hence the one-loop determination of the Landau pole can be misleading. A typical example 
is given by the two minimal DM cases in Eqs. (17)-(18). In this respect, we also pointed out 
the existence of a previously overlooked d = 5 operator which is responsible for a fast decay of 
(1, 7, 0)5, thus ruling out the scalar minimal DM candidate [44]. 

The final set of states which satisfy all the above constraints is collected in Tables 1-2. For 
these, we studied current bounds on their masses coming from their potential effects on BBN as 
well as from their production and detection at colliders. In particular, we found that for most of 
the states decaying through d = 5 operators and being thermally produced in the early Universe, 
their abundances are sufficiently diluted not to affect standard BBN. The notable exceptions 
are those uncolored cases, where the decay rates are either loop suppressed or proceed through 
long cascades leading to high-multiplicity final states (they are listed in Table 7). At colliders, 
the lightest particle of the multiplets in Tables 1-2 are, barring few exceptions, stable on the 
detector scale. For the color singlets we found that the current mass bounds are of few hundred 
GeV if the lightest particle of the multiplet is charged, whereas for neutral states the detection 
is more difficult and hence the bounds he below 100 GeV (cf. Table 8). On the other hand. 
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for colored multiplets the mass bounds strongly depend on the hadronization process and the 
nuclear interactions with the detector material. 

A this point, a natural question to ask is the following: What is accidental matter good 
for? Who ordered that? Besides the case of minimal DM, we note that the scalar multiplets in 
Tables 1-2 could easily improve the stability of the renormalizable Higgs potential.^® Close to 
the EFT cut-off, the potential can again be destabilized in cases where d = 5 operators exist 
containing only scalars, namely (1,5, 0)5, (1,5,1)5, (1,5, 2)5, (1,7, 0)5, (8,1,0)5, (8,1,1)5, 
(8, 3, 0)5 and (8, 3,1)5. However, at such large held values the whole tower of operators should 
be considered and stabilization is expected to be recovered via d > 6 operators. Finally, a 
charged thermal relic with ~ (10^ — 10^) s and abundance just below the D/H bound (which 
can be achieved for some of the accidental matter states listed in Table 7) may also help to 
resolve the standard BBN Lithium problems [74-79]. More generally, accidental matter should 
be seen as a purely phenomenological possibility. New physics might manifest itself in a way 
we were not expecting and thus the direct search strategies should cover diverse scenarios. 
In particular, the typical signature of accidental matter is the presence of charged/colored 
particles which are stable on the scale of particle detectors and which have no chances to be 
detected through indirect searches. Consequently, high-energy colliders will be the only means 
of probing such scenarios. New experiments in the near future (LHC-H, MoEDAL, etc) will 
have the capabilities to further explore their parameter space. 

We end by noting that an improvement in p-decay bounds by an order of magnitude [141] 
or failure to observe neutrinoless double beta decay with inverse neutrino mass hierarchy in the 
next generation of experiments [142] would put some pressure on this setup. In particular, (i) 
if neutrino oscillation experiments were to conhrm the inverse neutrino mass hierarchy, then 
our EFT setup predicts an observable neutrinoless double beta decay signal. Failure to observe 
one in the next generation of experiments would imply the presence of NP degrees of freedom 
below the EFT cut-off which couple to SM fermions. Finally, (ii) there is already a mild tension 
between the proton decay bounds and neutrino mass measurements, if one assumes a common 
EFT scale for both phenomena. This tension would be strengthened by future improvements 
in p-decay bounds or by a positive indication of a quasi-degenerate light neutrino spectrum, 
requiring a signihcant scale separation between the relevant L and B violating operators. 
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A Two-loop Landau Poles 


In this Appendix we provide the RG evolution of the gauge couplings and study the emergence 
of the associated Landau poles. In this way one can set an upper bound on the dimensionality 
of the extra representation, by requiring that no Landau poles are generated below Agg ~ 10^^ 
GeV (cf. the discussion in Sect. 2.4). 

The two-loop RG equation for the three gauge couplings gt {i = 1, 2,3), read 


d bij 
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where cni = ^ and t = ^ log The one- and two-loop beta function are [143] (no summation 
over i) 
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where Gi denotes the Gth gauge factor, S 2 and C 2 are the index (including multiplicity factors) 
and the quadratic Gasimir of a given (fermionic {F) or scalar (S')) irreducible representation; 
K = 1,1 for Dirac and Weyl fermions and p = 1, | for complex and real scalar helds, respectively. 
The Yukawa contribution in the two-loop beta function is neglected. In fact, the extra states we 
want to introduce do not couple with SM fermions, so that the Yukawa contribution does not 
grow with the dimensionality of the extra representation. Employing the GUT normalization 
for the abelian factor, we use the values ai{mz) = 0.016923, a 2 {mz) = 0.03374, and a^{mz) = 
0.1173 for the onset of the RG running [17, 144]. For simplicity, the extra state y is integrated 
in at mz = 91.188 GeV [17]. The scaling of the Landau pole with m-y- is approximately linear. 

For the cases where the SM is extended with a representation charged under SU(3)c and/or 
SU(2)x,, the results are summarized in Table 9, which provide a useful reference for estimating 
the bound on the dimensionality of the extra representations, by requiring that no Landau 
poles are generated below a given scale. 

The analysis has been repeated for all the states considered in this work, which can si¬ 
multaneously transform under SU(3)c and SU(2)y, and have a non-zero hypercharge as well. 
These include extra representations interacting with SM helds via d = 5 operators (but which 
cannot decay into SM states via renormalizable interactions), for which the results are reported 
in Table 10. Moreover, we investigated the Landau pole constraints for those extra scalar 
representations that can couple to the Higgs boson at the renormalizable level (cf. Table 5). 
Among them the only ones that do not couple to SM fermions at the renormalizable level and 
that survive the perturbativity criteria are the renormalizable cases of Table 1. In this respect, 
we mention a marginal case: (1,6, 1 / 2)5 for which = 6.6 x 10^^ GeV. Finally, we also 

checked the possibility of having multiplets decaying via d > 5 operators and whose neutral 
LP might be compatible with cosmological constraints. No cases beyond those of minimal DM 
(cf. Table 1) and with a Landau pole above 10^^ GeV are found. 
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RS 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

> mpi 

3> mpi 

3> mpi 

S> mpi 

S> mpi 

3> mpi 

1.4 X 10®® 

4.0 X 10® 

4.7 X 10® 

1.7 X 10^ 

8 

3> mpi 

3> mpi 

3> mpi 

8.3 X 10^2 

5.0 X 10® 

5.2 X 10® 

< 10® 

< 10® 

< 10® 

< 10® 

27 

1.3 X 10^ 

1.4 X 10^ 

< 10=^ 

< 10^ 

< 10® 

< 10® 

< 10® 

< 10® 

A 

o 

CO 

< 10® 

GS 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

3> mpi 

3> mpi 

3> mpi 

S> mpi 

3> mpi 

7.2 X 10^® 

2.0 X 10^ 

7.0 X 10^ 

4.7 X 10® 

1.1 X 10® 

3 

3> mpi 

3> mpi 

3> mpi 

> mpi 

1.1 X 10® 

4.7 X 10^ 

2.4 X 10® 

< 10® 

A 

o 

CO 

A 

o 

CO 

6 

3> mpi 

3> mpi 

3> mpi 

5.5 X 10^ 

1.6 X 10^ 

1.1 X 10® 

< 10® 

< 10® 

A 

o 

CO 

A 

O 

CO 

8 

3> mpi 

3> mpi 

> mpi 

1.1 X 10® 

4.2 X 10® 

< 10® 

< 10® 

< 10® 

A 

o 

CO 

A 

O 

CO 

10 

3> mpi 

3.8 X 10^ 

8.5 X 10^ 

1.3 X 10^ 

< 10® 

A 

o 

CO 

< 10® 

< 10® 

< 10® 

A 

O 

CO 

15 

3> mpi 

7.3 X 10^ 

2.0 X 10^ 

< 10® 

< 10® 

< 10® 

< 10® 

< 10® 

< 10® 

A 

O 

CO 

15' 

1.3 X 10^ 

< 10^ 

< 10^ 

< 10® 

< 10® 

< 10® 

< 10® 

< 10® 

A 

o 

CO 

A 

O 

CO 

24 

1.9 X 10^ 

< 10^ 

< 10^ 

< 10® 

< 10® 

< 10® 

< 10® 

< 10® 

< 10® 

< 10® 

WF 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

> mpi 

3> mpi 

3> mpi 

S> mpi 

8.3 X 10®^ 

5.0 X 10® 

3.7 X 10® 

1.3 X 10^ 

2.3 X 10® 

A 

o 

CO 

8 

3> mpi 

3> mpi 

7.3 X 10^ 

3.2 X 10^ 

1.4 X 10® 

< 10® 

< 10® 

< 10® 

< 10® 

< 10® 

27 

< 10^ 

< lO^ 

< 10=^ 

< 10® 

< 10® 

< 10® 

< 10® 

< 10® 

< 10® 

< 10® 

DF 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

3> mpi 

3> mpi 

3> mpi 

8.5 X 10^® 

9.4 X 10^ 

8.4 X 10^ 

4.4 X 10® 

1.0 X 10® 

A 

o 

CO 

< 10® 

3 

3> mpi 

3> mpi 

8.1 X 10^^ 

5.4 X 10® 

4.3 X 10® 

< 10® 

< 10® 

< 10® 

< 10® 

A 

o 

CO 

6 

3> mpi 

1.9 X 10^^ 

6.4 X 10^ 

3.3 X 10® 

< 10® 

A 

o 

CO 

< 10® 

< 10® 

A 

o 

CO 

A 

o 

CO 

8 

3> mpi 

1.5 X 10^ 

1.0 X 10^ 

1.3 X 10® 

< 10® 

< 10® 

< 10® 

< 10® 

A 

o 

CO 

A 

o 

CO 

10 

2.6 X 10^ 

< 10^ 

< 10^ 

< 10® 

< 10® 

A 

o 

CO 

< 10® 

< 10® 

< 10® 

< 10® 

15 

3.2 X 10^ 

A 

o 

CO 

< 10^ 

< 10® 

< 10® 

< 10® 

< 10® 

< 10® 

A 

o 

CO 

A 

o 

CO 


Table 9: Two-loop Landau poles (GeV) for the SM augmented with an extra multiplet (integrated in at mz) which is charged 
under SU(3)c and/or SU( 2 )i {Y = 0). The rows and columns denote respectively the SU(3)c and SU( 2 )l representations, while 
the four subtables correspond to the cases of an extra real scalar (RS), complex scalar (CS), Weyl fermion (WF) and Dirac fermion 
(DF). mpi ~ 10^® GeV is the Planck mass. The cases where the two-loop Landau pole is below 10^® GeV are emphasized in red. 







































Spin 

X 

AL^dllGeV] 

0 

(1,2, 3/2) 

> mpi (gi) 

0 

(1,2, 5/2) 

> mpi {gi) 

0 

(1,5,0) 

> mpi (gi) 

0 

(1,5,1) 

> mpi (gi) 

0 

(1,5,2) 

3.5 X 10^® (^i) 

0 

(1,7,0) 

1.4 X 10^® (c/a) 

0 

(3,1,5/3) 

> mpi {gi) 

0 

(3, 2, 5/6) 

> mpi {gi) 

0 

(3,2,11/6) 

5.5 X 10^9 (c/i) 

0 

(3, 3, 2/3) 

> mpi (gi) 

0 

(3, 3, 5/3) 

3.2 X 10^" (c7i) 

0 

(3,4,1/6) 

> mpi {g2) 

0 

(3,4, 5/6) 

> mpi {g2) 

0 

(6, 2,1/6) 

> rnpi {gi) 

0 

(6, 2, 5/6) 

> rnpi {gi) 

0 

(6, 2, 7/6) 

> mpi {gi) 

0 

(6,2,11/6) 

4.0 X 10^2 

0 

(6,4,1/6) 

5.5 X 10^ (c/a) 

0 

(6,4, 5/6) 

5.0 X 10^ (c/a) 

0 

(8,1,0) 

> mpi {gi) 

0 

(8,1,1) 

> mpi {gi) 

0 

(8,3,0) 

> mpi {gi) 

0 

(8,3,1) 

1.0 X lO^" (c7i) 

0 

(27,1,0) 

1.3 X 10^ (c73) 

1/2 

(1,3,2) 

1.4 X 10^3 

1/2 

(1.4.1/2) 

8.1 X 10^® (^a) 

1/2 

(1,4, 3/2) 

2.7 X 10^® (c7i) 

1/2 

(3, 3,4/3) 

9.3 X 10^° (c/i) 

1/2 

(3, 3, 5/3) 

1.6 X 10® (c/i) 

1/2 

(3,4,1/6) 

5.4 X 10® (c/a) 

1/2 

(3,4, 5/6) 

5.3 X 10® (c/a) 

1/2 

(3,4, 7/6) 

5.2 X 10® (^a) 

1/2 

(6,1,1/3) 

> (di) 

1/2 

(6,1,2/3) 

> (di) 

1/2 

(6, 2,1/6) 

1.9 X 1012 (^3) 

1/2 

(8,1,1) 

4.0 X lOi® (^1) 

1/2 

(8, 2,1/2) 

1.5 X lOi' (C73) 

1/2 

(15,1,1/3) 

3.2 X 10® ((73) 

1/2 

(15,1,2/3) 

3.2 X 10® ((73) 

1/2 

(15,2,1/6) 

3.6 X 102 (^3) 


Table 10: List of extra multiplets which can decay into SM particles via d = 5 operators 
(states decaying via d = 4 operators have been already subtracted) and corresponding two-loop 
Landau poles evaluated by integrating in the new states at niz (helds with zero hypercharge 
are understood to be real). In the third column, the symbol in the bracket stands for the gauge 
coupling responsible for the emergence of the Landau pole. 
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B S\J{2)l decompositions 


By denoting the generators in the fundamental representation of SU( 2 ) 2 , as T“ = ct“/ 2 (with 
being the Pauli matrices and a = 1, 2, 3), we dehne their action on the (2j + l)-dimensional 
completely symmetric tensor Xiii 2 -i 2 j (*i) * 2 , • • •, = 1 , 2 ) as 

^ ixiih—hj) ~ '^iik Xki2...i2j + Xiik...i2j + • • • + Xiih-.-k ■ (38) 

In general, we arrive at the following embedding of the properly normalized eigenstates: 


X11...1 

X11...2 


1 

a / B2j,0 




X12...2 

X22...2 


^jB2j,2j-l 
^ 1 y-t 


-i+1 




(39) 


where the superscripts denote the eigenvalue, is the binomial factor ^ and 

the normalization of the states is such that 


= ivr+ly-'r+...+ \x-’*x+\x-^ 


Let us consider, for instance, the case of the 811(2)^ Higgs doublet: 


H,=H^ 

H2=Ho, 


(40) 


(41) 


where the electric charge eigenstates are obtained through the formula Q = + Y. In par¬ 
ticular, since in the unitary gauge: = 0, Imhfo = 0 and Rehfo = :^('^ + ^)) whenever 

the Higgs doublet is contained in the effective operator responsible for the decay of y, it might 
happen that not all of the components of y can directly decay through the effective operator. In 
the following, we provide the SU(2)i decomposition for the three uncoloured multiplets whose 
decay, depending on the mass spectrum, might proceed via off-shell cascades (cf. Table 6 ): 


y = ( 1 , 2 , 5 / 2)5 


• SU( 2)5 embedding: 


Xl — X+3 

X 2 = y +2 • 

(42) 

• Operator: 


0 , = . 

(43) 

• Decomposition in the unitary gauge: 


Oi = ^ xUe^'e^iv + h). 

(44) 


Notice that y +3 does not couple directly to SM particles. Hence, whenever it is the LP it 
will decay through an off-shell emission of x +2 (cf- Fig- 6 ). 
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X= ( 1 , 5 , 1)5 

• SU(2)i;, embedding: 


Xiiii — X +3 

X1112 = 

X1122 = :j^X+i 
X1222 = ^Xo 
X2222 = X-i ■ 


(45) 


• Operator: 

O, = . (46) 

• Decomposition in the unitary gauge: 

0,^\xl(.v + hY. (47) 

Notice that only xo can directly decay through Oi. If xo is not the LP in the multiplet, 
the charged LP will cascade decay through off-shell components which end up into xo- 


X= (1,5, 2)5 

• SU(2)i embedding: 


Xiiii — X +4 
X1112 = ;^X +3 

X1122 = ;^X+2 ( 48 ) 

X1222 = ;^X+i 
X2222 = Xo • 


• Operator: 

Oi = . (49) 

• Decomposition in the unitary gauge: 

Oi = ^xl{v + hY. (50) 

Similarly to the previous case, only xo can decay through Oi, while the charged compo¬ 
nents decay through off-shell cascades. 
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